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ABSTRACT 

In  order  to  establish  a  model  which  would  more  accurately  predict 
discharge  characteristics  of  a  cell  at  any  arbitrary  current,  and  which 
would  better  fit  the  family  of  discharge  curves ,  testing  was  done  based 
on  the  Shepherd  equation  using  a  transparent  lead-acid  battery. 

A  set  of  discharge  curves  at  .6,  1.5,  3-6,  and  5-^+A,  preceded  by  a 
3.6A  preconditioning  discharge  to  remove  the  "memory  effect,"  was  obtained. 
The  battery  was  cycled  continuously  through  deep  discharges  in  an  attempt 
to  determine  secular  changes  during  the  life  of  the  battery. 

Several  modifications  to  the  Shepherd  equation  were  proposed  in 
order  to  obtain  a  better  fit  to  discharge  data.   These  included  modeling 
the  capacity  as  a  function  of  current  as  Q=CI^-n',  Peukert's  equation; 
modeling  the  ohmic  polarization  as  a  linear  function  of  accumulated  charge; 
and  modeling  the  diffusion  polarization  as  a  function  of  accumulated  charge 
alone  and  not  current . 

Conclusions  were  drawn  on  the  suitability  of  Shepherd's  equation. 
The  inclusion  of  the  changed  ohmic  polarization  term  gave  a  good  fit  for 
individual  curves  which  also  suggested  the  change  in  the  diffusion  polari- 
zation term.   Since  a  small  current  density  approximation  was  used  in 
Shepherd's  equation,  developing  a  model  based  on  the  original  terms  was 
recommended.  Also  recommended  were  the  modification  of  a  model  to  predict 
characteristics  during  a  cycle  with  various  discharge  currents  and  further 
life  cycle  testing  to  determine  the  secular  changes  in  the  battery. 
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CHAPTER  1 
INTRODUCTION 
1 . 1   BACKGROUND 

Predicting  the  transient  and  long-term  characteristics 
of  a  cell  under  arbitrary  charging  and  discharging  histories 
is  a  problem  confronting  a  system  designer  which  could  be 
overcome  if  there  were  an  electrochemical  model  of  the 
lead-acid  cell.   With  this  model,  it  would  theoretically 
become  possible  to  determine  the  optimum  scheme  for  getting 
energy  into  and  out  of  a  battery. 

To  what  degree  is  the  efficiency  of  a  battery  affected 
by  different  types  of  discharge?   Should  the  designer  select 
a  pulse-type  discharge  over  a  steady  discharge?   What  effect 
does  this  have  on  the  system  as  a  whole?   To  what  degree  does 
the  age  of  the  battery  affect  the  efficiency?   Low  current  or 
high?   What  exactly  is  the  effect  of  temperature  and  current 
on  the  charging  and  discharging  history  of  the  arbitrary 
battery?  Will  current  models  handle  the  problem  of  the  age 
of  the  battery?   Does  the  model  provide  a  fit  merely  to  a 
single  discharge  curve  or  for  the  family  of  curves?   Does  it 
account  for  the  condition  of  the  battery?  What  effect  has 
gassing  on  the  capacity  for  full  charge  or  discharge? 
Questions  such  as  these  pose  a  great  problem  for  the  designer 
of  a  model. 

Naturally,  there  are  numerous  other  criteria  which  might 
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affect  the  optimum  charge/discharge  history  of  a  battery,  so 
it  becomes  very  important  to  select  a  battery  type  to  test 
which  possesses  characteristics  which  are  relatively  similar 
to  other  cells. 

The  lead-acid  cell  was  chosen  for  this  investigation  for 
several  reasons.   Lead-acid  batteries  are  the  most  widely 
used  batteries  in  the  world,  and  for  that  reason,  they 
represent  the  greatest  commercial  interest.   Portable  lead- 
acid  batteries  have  automotive  use,  use  in  motorcycles, 
railroads,  all  types  of  military  vehicles  on  the  land,  sea, 
and  air,  and  additional  uses  not  related  to  vehicles. 

Additionally,  there  is  a  wealth  of  specific  information 
concerning  lead-acid  batteries  in  literature.   Volume  after 
volume  has  been  written  exploring  the  numerous  characteristics 
of  this  type  of  battery  in  great  detail.   Chapters  upon 
chapters  have  been  written  about  the  history,  electrochemistry, 
applications  and  performance  characteristics  of  lead-acid 
batteries.   Detailed  studies  have  been  done  on  performance  and 
structural  aspects,  as  well  as  on  self-discharge  effects. 

Because  the  author  chose  as  the  aim  of  his  thesis,  the 
establishment  of  an  electrochemical  model  of  the  cell,  based 
on  Shepherd's  equation  (1965),  it  was  necessary  to  chose  a 
battery  type  for  the  experimentation  which  was  rechargable, 
lightweight,  and  portable,  in  addition  to  having  similar 
characteristics  of  other  types  of  cells.   The  lead-acid 
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battery,  in  addition  to  being  rechargable,  is  available  in  a 
clear  plastic  case,  thus  affording  the  investigator  a  visual 
means  to  determine  some  battery  conditions  which  are  not 
readily  available  empirically  in  a  battery  with  an  opaque  case. 

Also,  it  was  felt  that  the  characteristics  of  some  of 
the  other  types  of  cells  are  similar  to  those  of  the  lead-acid 
cell  and  that  a  model  for  the  one  may  be  applicable  to  the 
others. 

1.2   HISTORY  OF  LEAD-ACID  BATTERIES 

The  original  incentive  toward  the  development  of  an 
electrochemical  source  of  current  was  the  "galvanic  cell" 
which  was  based  on  Galvani ' s  "frog  leg"  experiment,  circa 
1789. 

Volta,  after  1792,  developed  several  useful  current 
sources,  such  as  the  Volta  pile.  He  also  formulated  the 
first  electrochemical  theory  on  the  subject. 

Ritter  developed  a  cell  in  about  1803,  Ritter's  pile, 
which,  when  disconnected  from  an  external  source  of  current, 
would  yield  a  current  in  the  opposite  direction.   This  was  the 
basis  for  the  secondary  cell  or  storage  battery. 

In  1860,  Plante  demonstrated  his  cell  in  which  two  lead 
sheets  separated  by  rubber  were  rolled  into  a  spiral  and 
placed  in  a  10%  sulfuric  acid  solution.   When  charged,  the 
amount  of  energy  depended  upon  the  amount  of  lead  oxide  formed. 
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By  1880,  W.V.  Siemens  had  invented  the  dynamo,  which 
was  then  installed  in  central  power  stations .   Batteries  were 
used  for  load  leveling  purposes,  being  charged  with  energy 
produced  at  night. 

Around  18  81,  Faure  coated  the  lead  foil  with  lead  oxide 
for  the  purpose  of  accelerating  Plante's  forming  process. 
He  used  various  lead  current  collector  grids. 

1883  brought  us  the  Plante  plate  which  was  actually 
designed  by  Tudor.   In  this  design,  he  pasted  lead  oxide  onto 
cast  ribbed  pure  lead  plates.   Grid  plates,  plates  cast  from 
antimonial  alloys  were  also  developed  at  about  the  same  time. 

From  1890  to  1920,  wood  was  used  to  separate  the  plates 
and  then  rubber.   Around  1927,  porous  ebonite  and  microporous 
rubber  was  used  as  separators. 

Around  1965,  batteries  were  designed  to  be  relatively 
maintenance-free  for  use  in  portable  items,  and,  with  the 
advent  of  plastic  battery  cases  and  lightweight  construction, 
automation  and  mass-production  became  feasible. 

The  1970s  gave  us  new  applications  of  the  lead-acid 
battery,  with  increasing  interest  in  Traction  batteries  coming 
to  the  forefront.   Water-activated  and  maintenance-free 
automobile  batteries  were  designed.      Kordesch,  (.1977) 

These  historical  developments,  some  important,  some 
merely  historically  interesting,  represent  the  basis  for  the 
lead-acid  battery  as  it  is  known  today. 
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1 . 3   THE  NEED  FOR  A  MODEL 

A  system  designer  needs  an  accurate  electrochemical 
battery  model  in  order  to  produce  a  design  which  makes  optimum 
use  of  these  batteries .   With  the  use  of  computer-aided 
design,  a  functional  description  of  various  battery  character- 
istics would  greatly  simplify  the  present  cumbersome  task  of 
establishing  an  algorithm  for  interpolation  of  battery 
characteristics  from  its  family  of  discharge  curves.   An 
algorithm  for  predicting  the  response  to  an  arbitrary  (non- 
constant  current)  discharge  or  charge  would  be,  at  the  very 
least,  extremely  difficult,  if  not  prone  to  significant 
inaccuracies. 

Several  empirical  models  have  been  developed  in  the  past 
to  show  the  dependence  of  capacity  on  discharge  current  I 
and  discharge  time  t  {Schroder  (1894)  ,  Liebenow  (1897)  , 
Peukert  (1897),  Rabl  (1936),  Davytan  (1946),  and  Shepherd 
(1965)}.   Shepherd's  (1965)  model  (discussed  in  detail  in 
Chapter  II)  gives  a  relationship  between  cell  voltage  on  dis- 
charge as  a  function  of  current  density  and  time  based  on  the 
theoretical  laws  of  electrochemistry.   Shepherd's  equation 
(1965)  serves  as  the  basis  for  this  investigation  into  estab- 
lishing a  State  Model. 

Very  early  in  this  investigation,  it  was  found  that  the 
equation  provided  a  satisfactory  fit  to  a  single  discharge 
curve  (complete  discharge  at  constant  current)  but  a  poor  fit 
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for  the  family  of  discharge  curves.   A  typical  family  of 
discharge  curves  is  shown  in  figure  1-1.   Also,  Shepherd 
did  not  include  any  term  which  would  account  for  secular  life 
changes  in  the  cell. 
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Fig.  1-1.   Discharge  Curves  for  a  Typical  100-Ah  (20-h  rate) 
SLI  Lead-Acid  Battery. 


The  emphasis  of  this  thesis  was  to  establish  a  model 
which  would  more  accurately  predict  the  discharge  characteris- 
tics of  a  cell  at  any  arbitrary  current,  thereby  providing  a 
better  fit  for  the  family  of  discharge  curves.   Basing  his 
investigation  into  a  State  Model  on  the  Shepherd  equation 
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(1965) ,  the  author  has  chosen  to  call  his  model  the  "Modified 
Shepherd  Equation."   It  was  also  the  desire  of  the  author 
to  experimentally  determine  the  secular  life  changes  of  the 
lead-acid  cell  and  to  include  a  model  of  these  changes  in  the 
Modified  Shepherd  Equation. 

The  problem  of  needing  an  adequate  model  was  made  clear 
to  the  author  while  working  on  a  Diesel  submarine  computer 
synthesis  model  at  the  Naval  Ships  Engineering  Center  in 
Washington,  D.C.   A  great  deal  of  time  was  invested  in  estab- 
lishing a  model  of  the  lead-acid  cell  which  would  adequately 
predict  the  capacity  of  the  cell  at  arbitrary  discharge  cur- 
rents.  This  was  accomplished  simply  as  a  curve  fit  to  the 
capacities  determined  from  the  family  of  discharge  curves. 

Due  to  the  complexity  of  the  problem,  no  attempt  v/as  made 
to  use  the  discharge  characteristics  in  that  model.   If  a  model 
such  as  the  one  proposed  in  this  thesis  had  been  available,  not 
only  would  a  significant  number  of  man-hours  have  been  saved, 
but  also  another  degree  of  flexibility  would  have  been  allowed 
in  the  synthesis  model. 

Also,  the  problem  of  the  design  of  an  electric  vehicle  is 
essentially  identical  in  this  respect.   Moreover,  the  battery 
manufacturer,  by  supplying  only  a  few  coefficients  for  a 
battery  model  would  be  providing  the  characteristics  of  his 
product  in  a  very  concise  form  which  would  require  no  inter- 
polation for  the  user's  specific  application. 
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1.4   BATTERY  USAGE 

Since  the  earliest  "galvanic  cell",  the  lead-acid  storage 
battery  has  been  in  demand  for  numerous  purposes.   Marine 
uses  alone  are  diverse. 

Marine  batteries,  which  are  rugged  and  durable,  have 
given  long  service  to  ships  of  the  Navy  and  Merchant  Marine. 
These  batteries  function  under  extremely  adverse  conditions 
and  are  subject  to  continual  vibration. 

These  batteries  have  been  used  in  all  types  of  ships  from 
the  largest  aircraft  carrier  to  the  smallest  pleasure  craft  to 
power  radios,  lights,  pumps,  hoists,  and  compressors.   They 
are  used  in  torpedoes  and  to  fire  guns.   They  are  used  in  the 
kitchens  of  ships  to  provide  power  to  stoves  and  refrigerators. 

Batteries  power  alarms  and  emergency  lighting  and  even 
pump  bilges  and  hoist  anchors.   Many  are  used  in  ignition 
circuits  of  internal-combustion  engines  of  small  craft.   They 
provide  emergency  power  for  submerged  Nuclear  submarines  and 
the  propulsion  for  submerged  Diesel  submarines. 

Marine  batteries  are  used  for  gyrocompass  emergency  power 
and  for  telephones,  and  emergency  generator  start-up/  or 
field  flashing. 

For  many  years,  scientists  and  engineers  have  studied  the 
possibility  of  converting  from  gasoline  powered  motor  vehicles 
to  electrically  powered  motor  vehicles,  from  hybrids  to  all- 
electric  delivery  vans.   Much  interest  was  in  all-electric 
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automobiles.   It  has  been  found  that  there  are  no  savings  in 
energy  if  batteries  simply  replace  the  gasoline  engine. 
Electrically  powered  vehicles  would  offer  a  significant 
pollution  advantage  over  gasoline  and  Diesel  engines. 
Batteries  could  be  used  to  provide  additional  energy  require- 
ments in  hybrid  gasoline/electric  vehicles  for  load  leveling 
purposes  or  for  applications  where  the  possible  recovery  of 
energy  through  regenerative  breaking  is  significant. 
(Kordesch  1977) 

Batteries  have  been  classified  according  to  usage, 
according  to  Bode  (1977)  : 

1.  SLI  Batteries  (Starting,  Lighting  and  Ignition) 
These  batteries  are  used  to  start  an  internal- 
combustion  engine  and  to  provide  power  to  the 
electrical  system  when  the  engine  is  not  running. 
The  largest  number  of  storage  batteries  are  for 

use  in  the  automobile  and  are  used  for  this  purpose 
more  than  for  any  other. 

2.  Stationary  Batteries 

These  generally  provide  emergency  power  when 
the  primary  source  of  power  is  not  available. 
These  include  such  applications  as  communications 
systems,  electric  utilities,  computer  systems, 
emergency  lighting,  and  railways  to  provide  peak 
loads,  emergency  power  or  smoothing  of  ac  power. 
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3.  Motive  Power  (Traction)  Batteries 

Lead  storage  batteries  provide  power  for  the 
propulsion  of  devices  such  as  electric  lift  trucks, 
mining  equipment,  street  delivery  vehicles,  and 
other  types  of  material  handling. 

4.  Special  Purpose  Batteries 

Some  examples  of  batteries  in  this  category 
are  aircraft,  submarine,  special  military  and  small 
sealed  batteries  for  consumer  applications.   Some  of 
these  uses  have  already  been  discussed. 
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1.5   BATTERY  LIFE 

The  life  of  a  battery  can  be  defined  as  the  period  of 
time  during  which  it  has  the  ability  to  perform  within  certain 
specified  limits.   In  other  words,  the  end  of  life  occurs 
when  the  battery  can  no  longer  perform  within  these  limits. 
These  limits  actually  depend  upon  the  application  and  the  re- 
quirements of  the  system  they  serve.   The  end  of  life  of  an 
automobile  battery  occurs  when  it  no  longer  is  able  to  provide 
sufficient  cranking  power  to  start  the  engine.   For  a  battery 
powered  vehicle,  it  may  be  when  the  vehicle  will  no  longer 
travel  a  certain  number  of  miles  at  a  specified  speed;  in 
traction  batteries,  the  decrease  of  capacity  below  a  certain 
value  may  be  the  criterion  on  which  the  end  of  life  decision  is 
based.   Generally,  however,  systems  are  designed  to  meet  all 
specified  requirements  as  long  as  the  battery  is  able  to 
deliver  80%  or  more  of  its  new  capacity. 

Bode  (1977)  states  that  the  end  of  battery  life  is  best 
evaluated  in  the  laboratory  using  fixed  standards  and  speci- 
fied test  conditions.   The  ideal  requirement  is  that  the  life 
test  correspond  to  the  load  conditions  of  the  battery  in  the 
field,  however,  as  in  the  case  of  automobile  batteries,  ideal 
testing  in  the  laboratory  is  impossible  because  it  is  im- 
possible to  duplicate  such  diverse  conditions  as  automobile 
type,  driving  methods,  maintenance  and  environment. 
Figure  1-2  is  a  plot  of  the  change  in  capacity  during  the 
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lifetime  of  a  lead-acid  battery.   Very  little  information  is 
available  which  provides  details  on  the  change  of  characteris- 
tics of  a  cell  during  its  lifetime.   It  is  therefore  proposed 
that  this  data  be  obtained  by  experimental  means. 


Figure  1-2 
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Typical  Cell  Capacity  During  the  Lifetime  of 
a  Lead-Acid  Battery.   KORDESCH  (19  77) 


The  discharge  capacity  of  a  battery  has  been  found  to  be 
somewhat  dependent  upon  the  rate  of  previous  discharges.   Bode 
(1977)  calls  this  the  "memory  effect."   Vinal  (1955)  states 
that  "the  capacity  is  lower  if  the  discharge  has  been  preceded 
by  a  discharge  at  a  higher  rate,  and  it  is  higher  if  preceded 
by  a  discharge  at  a  lower  rate."   Therefore,  any  experimental 
procedures  should  take  the  "memory  effect"  into  consideration. 


b 
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1.6   REASONS  FOR  DECREASE  IN  LIFE  EXPECTANCY  OR  FAILURE 

According  to  Vinal  (1955) ,  storage  batteries  large  and 
small  and  in  widely  diverse  services  are  subject  to  conditions 
that  may  lead  to  excellent  service  or  through  accident  or 
abuse  to  something  less  than  satisfactory.   Under  his  heading 
"Sources  of  Trouble,"  he  writes  about  faulty  conditions 
which  if  detected  could  be  corrected.   Detection  and  correction 
of  faulty  conditions  might  prolong  the  satisfactory  battery 
life  into  years. 

He  also  states  that  as  a  result  of  studies  done  in  1947 
on  thousands  of  worn  out  batteries  which  had  been  used  in 
normal  service,  the  reasons  they  wear  out  became  well  known. 
In  comprehensive  data  relating  primarily  to  automobile 
batteries,  grid  corrosion  of  the  positive  plates  accounted  for 
42%  of  the  total  failures.   Naturally,  with  the  use  of  new 
alloys,  batteries  developed  since  that  time  should  have  less 
positive  plate  corrosion. 

Bode  (1977)  indicates  that  destruction  of  the  positive 
grids  occur  when  lead  sulfate  forms  a  heavy  layer  on  the  grid 
to  protect  it  from  being  dissolved  anodically.   One  rarely 
encounters  corrosion  of  the  negative  grid.   Hydrogen  corrosion 
occurs  in  the  presence  of  organic  substances  and  at  higher 
operating  temperatures.   The  result  of  corrosion  is  loss  of 
conductivity  and  grid  breakage. 
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Vinal  (1955)  indicated  that  7  to  10%  had  less  than  total 
life  span  due  to  shedding  of  positive  active  material, 
cracked  partitions  and  leaking  cases  accounted  for  nearly  one- 
third,  and  buckled  plates  and  short-circuited  separators  pro- 
duced failure  in  10  to  12%  of  the  batteries.   About  10%  of  the 
failures  were  caused  by  negative  plate  failure.   Sulfation 
which  causes  buckled  plates  and  negative  failures,  was  a  cause 
of  battery  failure  in  less  than  20%  of  the  batteries. 

Overcharging  loosens  active  material,  corrodes  positive 
grids,  causes  excessive  gassing  and  produces  sludging  of  a 
fine  brown  sediment.   It  may  cause  raised  temperatures  contri- 
buting to  battery  destruction  with  damage  to  plates  and  separ- 
ators, buckling  of  plates  and  needless  loss  of  water. 

Bode  (1977)  states  that  in  sulfation,  .  "hard"  lead 
sulfate  forms  if  chargeable  lead  sulfate  undergoes  long  or  high 
temperature  aging  process,  (recrystallization) .   Sulfation 
produces  loss  of  capacity  and  is  caused  by  long  standing  in  a 
discharge  condition,  too  high  acid  concentration,  increased 
self-discharge,  continuous  operation  between  45   and  50  C, 
and  prolonged  undercharging. 

He  also  indicates  that  sludging  of  the  positive  mass 
occurs  most  at  the  end  of  a  charge  and  at  the  beginning  of  a 
discharge.   Heavy  gassing  facilitates  sludging.   There  is  a 
greater  amount  of  sludging  in  discharges  with  large  current 
density.   The  type  of  discharge  influences  sludging.  If  acid 
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concentration  is  drastically  reduced,  temperature  increases  or 
current  density  is  diminished,  there  is  less  sludging  and 
the  battery  life  is  increased. 

Vinal  (1955)  concludes  that  occasional  overcharging  is 
beneficial,  but  habitual  overcharging  decreases  useful  battery 
life. 

He  indicates  that  undercharging  gradually  runs  down  cells 
producing  diminishing  values  of  specific  gravity.   Plates  be- 
come lighter  in  color  and  a  sludge,  a  fine  white  powder  which 
is  predominantly  lead  sulfate,  is  produced.   One  or  more  cells 
could  become  exhausted.   He  concludes  that  the  remedy  is  to 
charge  until  all  cells  are  in  normal  condition.   This  is  refer- 
red to  as  an  equalizing  charge.   Undercharging  is  the  most 
common  cause  of  buckling  of  plates  because  lead  sulfate  occu- 
pies more  space  than  the  original  material  producing  excessive 
strain  on  the  plates. 

Failure  of  the  charging  system  and  of  the  battery  ulti- 
mately is  blamed  on  the  battery.   When  an  automobile  reaches 
speeds  greater  than  or  equal  to  12  to  15  mph  the  battery  begins 
to  charge.   If  the  generator  doesn't  come  on  at  the  proper 
speed  the  battery  won't  receive  its  required  amount  of  charge. 
Also,  the  system  must  regulate  itself  properly  so  that  the 
battery  will  not  receive  massive  amounts  of  current  when  the 
automobile  reaches  highway  speeds.   Conditions  in  the  charging 
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system  which  drain  the  battery  and  cause  it  to  remain  dis- 
charged are  grounded  circuitry  and  partial  grounds. 
(Vinal,  1955) 

Other  possible  causes  of  shortened  battery  life  include 
corroded  terminals,  short  circuits,  separator  failures, 
worn-out  plates,  electrolyte  below  the  tops  of  the  plates, 
freezing,  when  the  battery  is  in  a  discharge  condition, 
impurities  in  the  solution,  spalling  of  the  active  material, 
growth  of  positive  plates,  shrinkage  of  negative  plates, 
leading,  and  battery  explosions.   The  interested  reader  is 
directed  to  the  classic  text  of  Vinal  (1955)  and  the  excellent 
texts  of  Bode  (1977)  and  Kordesch  (1977) . 
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CHAPTER  2 
THEORETICAL  ANALYSIS 
2.1   ELECTROCHEMISTRY  OF  THE  LEAD-ACID  BATTERY 

The  theory  of  operation  of  the  lead-acid  battery  is 
based  on  the  "double  sulfate  theory"  proposed  by  Gladstone 
and  Tribe,  (1881-1883)  . 

Pb02  +  Pb  +  2H2SO.»  ^1SCharg1  2PbS04  +  2H20        (2.1) 

charge 

The  principle  was  so  named  because  the  discharge  product 
of  both  the  lead  anode  electrochemical  oxidation  reaction  and 
the  lead  dioxide  cathode  electrochemical  reduction  reaction  is 
lead  sulfate. 

This  early  equation  reveals  the  stoichiometric  relation- 
ship of  reactants  and  products,  as  well  as  the  role  played  by 
the  sulfuric  acid. 

A  more  complete  representation  of  the  actual  electro- 
chemical reactions  is: 


ANODE:   Pb  +  HSOZ  dlschar9<!  PbSO„  +  H+  +  e         (2.2) 

charge 

CATHODE:   Pb02  +  3H+  +  HSO~  +  2e~— __— ±  PbSOi*  +  2H20     (2.3) 


OVERALL:   Pb  +  Pb02  +  2H+  +  2HSO4  »2PbS0i,  +  2H20     (2.4) 
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The  electrolytic  dissociation  of  sulfuric  acid  into 

H   and  HSO^  is  closer  to  the  true  picture  in  the  cell.   The 

+       —  2 
dissociation  into  2H   and  SOi*   occurs  in  solutions  much  more 

dilute  than  those  found  in  a  lead-acid  cell. 

During  overcharge,  the  electrolysis  of  water  occurs, 
yielding  hydrogen  and  exygen. 

Based  on  the  free-energy  change  of 

AG  =  -13.95  kcal/mole 
for  the  anodic  reaction 


(E   =  0.303  V  versus  H2/H+) 


and 


AG  =  -74.99  kcal/mole 


for  the  cathodic  reaction 


CE   =  -1.627  V  versus  H2/H+), 
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the   overall  reaction  values  are: 

AG  =  -88.94  kcal/mole 
(EQ  =  1.93  V) 

specific  energy  =  161  Wh/kg 

specific  capacity  =  83.5  Ah/kg 

Bode  (1977)  . 

In  practice,  significantly  lower  specific  energies  are 
exhibited,  but  the  theoretical  values  provide  an  upper  limit  to 
possible  future  developments.   For  a  complete  discussion  of 
electro-chemical  theory,  the  reader  is  directed  to  the 
excellent  text  by  Bode  (1977)  ,  and  the  classic  text  in  the  area 
by  Vinal  (1955) . 
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2.2   THE  SHEPHERD  EQUATION 

Shepherd  (1965)  has  proposed  a  static  model  for  cell 
discharge  characteristics  which  uses  simple  analytical 
expressions  of  chemicophysical  origin  and  describes  the  cell 
potential  E  as  follows: 

E(q,i)  -  Bo  -  VaCq)  -  V0  (i)  -  Vd(i,q)  (2_5) 

where  Eo=  open  circuit  cell  potential 

V  (q)  =  A  exp  (-  — £—   )  =  activation  polarization 
a  U 

Vq (i)  =  Roi  =  ohmic  polarization 

V, (i,a)  =  K  (  — -^ —  )  i  =  diffusion  (concentration) 
a         -  Q-q 


polarization 


i  =  the  discharge  current 

q  =  it  =  the  accumulated  charge 

Q  =  the  cell  capacity 

K  =  the  mean  polarization  resistance 

Ro  =  the  internal  cell  resistance 

Shepherd  used  "it"  in  place  of  "q"  in  his  equation 
and  Ro  corresponds  to  his  N.   The  activation  polarization 
term  is  the  single  term  which  does  not  have  any  particular 
theoretical  basis.   Its  inclusion  accounts  for  the  rapid 
drop  in  potential  at  the  beginning  of  discharge.   Figure  2-1 
shows  an  exagerated  view  of  this  phenomenon.   The  time 
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duration  of  this  initial  drop  was  found  to  be  in  order  of  se- 
conds.  Consequently,  the  description  of  cell  discharge  should 
not  be  significantly  affected  by  simply  ignoring  the  term;  this 
was  the  approach  used  by  Shepherd,  also. 

Another  phenomenon  shown  in  fig.  2-1,  not  treated  by  Shep- 
herd, is  known  as  the  "voltage  dip"  or  "coup  de  fouet"  where 
initially  the  voltage  drops  below  the  "steady  state"  value  and 
then  rises  to  this  level  within  a  few  minutes. 


Voltage  f 


"coup  de  fouet" 


Time 


Figure  2-1.   Cell  Voltage  versus  Time  for  a  Constant 
Current  Discharge 


The  ohmic  polarization  term  accounts  for  the  ohmic  resis- 
tance through  both  the  metal  and  the  solution.   This  has  been 
modeled  as  a  constant,  however,  it  is  felt  that  this  assump- 
tion may  be  erroneous. 
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The  specific  electrical  resistance  of  lead  oxide  is  on 
the  order  of  10" 6  Q/m   and  metallic  lead,  on  the  order  of 
10-7  fi/m.   The  electrolyte  has  a  range  of  specific  resistance 
of  about  125  fi/m  when  nearly  fully  charged  to  about  17  0  ft/m 
when  discharged.   The  product  of  the  two  reactions,  lead 
sulfate,  is  highly  resistive,  specific  resistance  on  the 
order  of  10 10  Q/m.      Based  on  these  values,  it  would  be  expect- 
ed that  the  ohmic  resistance  increase  as  the  reaction  proceeds 
and  becomes  fairly  significant  as  it  approaches  completion. 
The  diffusion  (concentration)  polarization  term  appears 
to  have  been  made  up  from  a  combination  of  what  Bode  (197  7) 
calls  the  charge  transfer  overvoltage,  n  ,  and  the  concen- 
tration overvoltage,  n  .   Following  Bode ' s  development, 
the  terminal  voltage  of  the  cell  is 


E   =   Erev  +  (nt  +  nc  }A  "  (rlt  +  Vc  "  IR        (2'6) 


E    is  the  open  circuit  voltage  and  IR  includes  the 
resistances  of  the  electrolyte  and  separator.   The  subscripts 
A  and  C  refer  to  the  anode  and  cathode  overvoltages ,  respec- 
tively.  All  of  the  overvoltages  decrease  the  cell  voltage  on 
discharge,  and  overvoltages  of  the  electrodes  are  related  to 
the  current  densities,  i.   Assuming  that  the  electrodes  have 
equal  surface  area  AA  =  A_,  =  A,  then  for  both  electrodes  the 
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following  expressions  apply: 


=  2-3  RT  log  i0  (1  -0  )  -  b  log  i  (2.7) 

zF 
and 

2.3RT  log  (1 i —  j  #  (2>8) 


n„  = 


zF  l, 


For  small  current  densities,  these  become 


tk  =  — (2.9) 


fc    zFiod  -0) 
and 


n  =  -5J i-  (2.10) 

zF     l. 


For  this  case,  one  of  the  terms  in  (2.6)  may  be  neglect- 
ed if  they  differ  by  an  order  of  magnitude.  The  terms  in  the 
equations  above  are: 

i0  =  the  apparent  exchange  density 

z  =  the  number  of  electrons  transferred  in  the 
determining  step 

F  =  the  Faraday  (measure  of  capacity) 

R  =  the  gas  constant 

T  =  the  temperature 

i   =  the  diffusion  limiting  current 
and 

0  =  the  degree  of  coverage.   (only  the  uncovered  portion 
(1-0)  of  the  electrode  can  react) 
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Replacing  the  constant  terms  in  n   by  a  single  constant  k 

yields 

n.  = (2.11) 

z       (1-0) 

If  the  uncovered  portion  of  the  electrode  is  related  to 
the  cell  capacity  Q  and  accumulated  charge  q  where  q  <  Q 
with  q  ■*■  Q  at  complete  discharge,  equation  (2.11)  becomes 


K     _  _     Q     .  (2.12] 

n,  -  -  jn.  

(1-0)         (Q-q) 


This  term  is  similar  to  Shepherd's  diffusion  polariza- 
tion but  with   the  striking  difference  that  it  is  not  a 
function  of  current.   If  one  were  to  fit  Shepherd's  equation 
to  a  single  discharge  curve,  this  difference  would  not  be 
apparent  because  the  product  of  K  and  i  in  equation  (2.5) 
would  be  equal  to  K  in  equation  (2.12) . 

Prior  to  the  discovery  of  the  above  relationship,  the 
author,  while  attempting  to  make  a  fit  for  a  family  of  dis- 
charge curves  to  Shepherd's  equation,  produced  fits  of  each 
curve  individually.   In  the  process  of  analyzing  these  results 
a  rather  surprising  discovery  was  made  in  that  the  product 
of  K  and  i  was  nearly  constant. 

Another  item  in  Shepherd's  equation  which  should  be 
addressed,  is  the  cell  capacity,  Q.   The  value  for  Q  should  be 
close  to  the  capacity  of  the  cell  at  the  particular  discharge 
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current.   In  fact,  it  needs  to  be  slightly  larger  than  the 
actual  capacity  to  avoid  predicting  infinite  overvoltages , 
i.e.   Q  >  q  for  all  values  of  q  at  the  particular  discharge 
current. 

The  most  widely  used  model  for  all  capacity  is  Peukert's 
(1897)  equation, 


Int  =  C  .  (2.13. 


In  this  equation  n  and  C  are  constants  which  can  be  determined 
by  using  the  discharge  rate  and  time  for  complete  discharge 
at  two  different  currents.   Since  the  cell  capacity  in  terms 
of  Q  is  the  desired  quantity,  the  substitution  of  t  =  Q/I  gives 


I(n-1)Q  =  C   or   Q  =  CI(l_n)  (2.14) 


The  values  of  n  and  C  can  be  calculated  as  follows: 


Qi  -  c  if1"  n) 
Q  =  C  I<1_n) 


(n-1) log  Ij  =  log  C  -  log  Qx 
(n-l)log  I2  =  log  C  -  log  Q2 


log  Q2  -  log  Qt 

n  -  +  1 

log  I j  -  log  I 2 
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The  value  for  n  is  in  the  range  of  from  1  to  2 .   For 
diminishingly  small  current  densities,  n  approaches  1  and 
for  high  current  densities,  n  approaches  2.   Typically,  the 
value  of  n  falls  v/ithin  about  1.3  to  1.4   This  equation  is 
not  based  on  any  particular  physical  laws,  but  serves  as  a 
good  empirical  fit  to  discharge  capacities. 

Baikie,  Gillibrand,  and  Peters  (1972)  expanded  on  this 
relationship  to  include  temperature  effects  in  the  range  of 
-18  C  to  +16°C;  for  current  densities  between  21  and 


165mA/cm2  it  has  the  form 


K0 (1  +  a6) 

Q  =  (2.15) 

z(n-i) 


Where  Q  is  the  effective  plate  capacity,  and  9   is  the  temp- 
erature in   C.   The  values  they  established  were 

n  =  1.4  for  both  plates 
K   =0.32  and  0.24  min, 
and  a   =   0.021  and  0.015  min/°C 


for  the  negative  and  positive  plates,  respectively. 

Even  if  these  equations  provide  a  perfect  correlation  to 
the  capacity  of  a  cell,  they  do  so  only  at  one  instant  and 
do  not  predict  changes  which  occur  throughout  the  life  of  the 
battery.   If  a  state  model  is  to  be  established  for  the  electro- 
chemical cell,  then  these  effects  must  also  be  included. 


39 
2.3    METHOD  OF  DETERMINING  COEFFICIENTS  FOR  THE  MODEL 

An  iterative  approach  using  a  digital  computer  v/as  used  to 
determine  the  coefficients  for  the  model.   This  curve  fitting 
method  basically  used  the  following  steps: 

1.  Choose  initial  values 

2.  Vary  each  value  independently 

n 
Trial  Value  =  Value  *  (Factor) 

where  n  =  -1,0,1   and  0   Factor   1 

3.  Choose  the  set  of  trial  values  which  gives  the 
smallest  difference  between  the  data  and  model. 

Steps  2  and  3  were  repeated  until  no  further  improvement 
was  required  by  the  user.   When  the  "best"  set  of  values  was 
found,  then  these  trial  value  were  used  in  the  next  iteration 
as  the  values.   When  a  particular  factor  failed  to  make  any 
improvement,  the  difference  between  the  factor  and  1  was  reduc- 
ed by  1/2  and  a  new  iteration  was  performed. 

For  runs  of  the  Shepherd  equation,  initial  values  were 
chosen  such  that: 


E  =  open  circuit  voltage 

Q  =  cell  capacity  plus  a  small  amount  to  prevent 
division  by  zero 


Ki  +  Ri  =  Es  -  discharge  voltage  at  t  =  0 
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This  method  finds  a  minimum  value  for  the  correlation, 
but  the  only  assurance  is  that  the  set  of  values  corresponds 
to  a  relative  minimum. 
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CHAPTER  3 
EXPERIMENTAL  PROCEDURE 
3.1   BATTERY  CHARACTERISTICS 

In  order  to  establish  a  model  which  would  more  accurately 
predict  the  discharge  characteristics  of  a  cell  at  any  arbi- 
trary current  and  throughout  its  life,  testing  was  performed 
using  a  lead-acid  battery. 

The  type  of  battery  used  was  a  low  capacity  motorcycle 
type  battery  in  a  clear  case.   The  low  capacity  was  used  so 
that  the  equipment  could  be  chosen  requiring  smaller  current 
carrying  capacities.   The  advantage  of  the  clear  case  over  an 
opaque  case  was  that  the  author  could  more  readily  observe 
conditions  of  the  battery,  including  the  amount  of  gassing 
during  charging,  the  level  of  the  electrolyte,  the  changes,  if 
any,  in  the  condition  of  the  plates  and  the  accumulation  of 
any  sludge,  or  fine-grained  sediment,  in  the  bottom  of  the  bat- 
tery case  without  destroying  the  battery  in  order  to  make  these 
observations . 

Table  3-1  lists  the  battery  specifications  and  overall 
dimensions.   Figure  3-1  shows  the  internal  configuration  of  a 
cell.   Each  cell  was  configured  with  anode  (-)  plates  on  the 
ends.   As  can  be  seen,  in  the  figure,  there  are  four  anode  (-) 
plates  and  three  cathode  (+)  plates  with  6  separators  in  the 
spaces  between  the  plates. 
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Figure  3-1.  Internal  Configuration  of  a  Cell. 


Manufactured  for: 


Battery  Type: 
JIS  Designation: 
Capacity: 
Case: 

External  Case  Dimensions: 


Wisco  Division 
ESB  Incorporated 
1222  18th  Ftreet 
Pacine,  WI  53403 

JA-11  6V       3  Cell 

6N6-3B 

6Ah  @  10-h  rate 

Clear  Plastic  Bottom 
with  Black  Plastic  Top 

L  =  10.7cm 
W  =  5.5cm 
H  =  11.0cm 


Table  3-1.  Battery  Specifications-  and  Dimensions. 


^8-> 

■mm. 
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Anode  Plate  #1 
Cell  #1  Battery  #1 


i»6 
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8  W7H< 


Connector 


Note:   All  plates 
appear  to  have  the  same 
dimensions . 


Area 


fU 


29cm 


Note:   This  tab  sup- 
ports the  plate  on  a 
rib  in  the  case  run- 
ning the  length  of  the 
cell . 


Figure  3-2.   Side  View  of  Anode  Plate  #1  as  Seen  from  the 
End  of  the  Battery. 


Table  3-2  lists  the  measurements  on  cell  #1  to  de- 
termine the  thicknesses  of  the  plates  and  separators.   These 
measurements  were  made  by  placing  a  rule  against  the  side  of 
the  case.   Some  degree  of  uncertainty  in  the  measurement  of  the 
plate  thicknesses  can  be  attributed  to  the  method  of  measurement 
required. 

Figure  3-2  shows  the  side  view  of  anode  plate  #1  which  is 
the  view  seen  when  looking  at  the  end  of  the  battery.   All 
plates  have  essentially  the  same  dimensions.   The  area  of  the 
plates  is  approximately  29cm2. 
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COMPONENT          DISTANCE   THICKNESS 

(mm) 

(mm) 

Anode  Plate  #1 

1.5 

1  .5 

Separator  #1 

4.0 

2.5 

Cathode  Plate  #1 

6.2 

2.2 

Separator  #2 

9.0 

2.8 

Anode  Plate  #2 

10.3 

1.3 

Separator  #3 

13.0 

2.7 

Cathode  Plate  #2 

15.0 

2.0 

a 

Separator  #4 

17.7 

2.7 

Anode  Plate  #3 

19.5 

1.8 

Separator  #5 

22.0 

2.5 

Cathode  Plate  #3 

24.0 

2.0 

Separator  #6 

26.5 

2.5 

Anode  Plate  #4 

28.2 

1.7 

AVERAGE  THICKNESS 


Anode 

1.575mm 

S.  Dev. 

0.222 

Var. 

0.0369 

Cathode 

2 . 07mm 

S.  Dev. 

0.115 

Var. 

0.0089 

Separators 

2 . 6  2mm 

S.  Dev. 

0.133 

Var. 

0.015 

Table  3-2.   Measurements  of  the  Distances  of  the  Components 
from  Edge  of  Anode  Plate  #1  and  the  Thicknesses 
of  the  Anode  and  Cathode  Plates  and  Separators. 
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hole 
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Figure  3-3.   Top  View  of  Test  Battery 


Table  3-3  shows  the  volume  of  the  materials  in  a  cell, 
including  the  deduced  volume  of  electrolyte  of  SOcc. 

Access  to  the  connectors  between  cells  #1  and  #2,  and 
cells  #3  and  #4  was  made  by  drilling  into  the  top  of  the  case. 
The  location  of  the  holes  is  shown  in  figure  3-3. 

The  following  observations  were  made  of  the  battery  in  its 
unused  condition.   The  anode  plate  viewed  from  the  side  had  a 
light  grey  appearance;  viewed  from  the  end  it  was  darker  grey. 
Bubbles  were  seen  on  the  plate  and  periodically  several  would 
rise  to  the  surface.   Cathode  plates  could  only  be  observed 
from  their  ends  due  to  the  configuration  of  the  cells.  They  had 
a  dark  brown  rusty  coating  on  their  surfaces .   Bubbles  were  not 
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INTERNAL 

VOLUME  (per  c 

ell)  : 

L  =  2  9mm 

W  =  51 . 5mm 

Height  to  upper  electrolyte 

level 

(marked  on  case 

)  = 

8  0mm 

Volume  - 

120cc 

VOLUME  OF  ELECTROLYTE: 

Vol.  of 

anode  (per  ce 

:11)  = 

2  9cm  * 

.1575cm 

*  4  = 

18. 

3cc 

Vol.  of 

cathode  (per 

cell) 

=  29cm 

*  .207cm  *  3  = 

18 

.  Occ 

Internal  volume  -  p] 

.ate  volume  = 

120cc  - 

36cc  = 

84 

cc 

Assume  solid 

volume 

;  of  separators 

=  4cc 

VOLUME  OF 

ELECTROLYTE  « 

=  80cc/ce 

ill 

Table  3-3.   Volume  of  Materials  within  a  Cell 
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generally  observed  on  these  plates. 

These  observations  were  made  several  weeks  after  the 
battery had  first  been  activated,  but  before  it  had  been 
used.   When  the  battery  had  first  been  activated,  the  anode 
and  cathode  plates  had  more  nearly  the  same  appearance,  which 
was  light  grey. 

The  separators  appeared  to  be  made  of  a  "fiber  wool" 
placed  on  a  paper  (or  felt)  backing.   The  backing  was  always 
placed  against  the  anode  plates,  while  the  "wool"  was  against 
the  cathode  plates.   Small  bubbles  were  entrapped  in  the  fiber. 

The  cavity  below  the  plates  had  a  very  slight  amount  of 
rusty-colored,  fine-grained  sediment. 

3.2   EXPERIMENTAL  SET-UP 

Figure  3-4  is  a  circuit  diagram  of  the  experimental  set-up 
using  a  double  pole-double  throw  switch.   The  circuit  was 
configured  in  this  manner  so  that  no  disconnecting  and 
reconnecting  would  be  required  between  charges  and  discharges 
and  a  single  ammeter  could  be  used  in  either  mode. 

Two  different  battery  charging  set-ups  were  used.   One  was 
simply  the  use  of  a  commercially  available  6  Volt  Tapering-Type 
battery  charger.   The  battery  charger  was  manufactured  by 
Mallory,  Indianapolis,  Indiana  (Engineering  Projects  Lab 
Number  1167),  Model  6SAC10. 

Two  chart  recorders  were  used.   One  was  a  4  channel  chart 
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Figure  3-4.   Circuit  Diagram  of  the  Experimental  Set-Up 
Using  a  Double  Pole-Double  Throw  Switch. 


recorder  manufactured  by  Gould   Brush  2400  Model  2007-4404-00 
Serial  00473.   All  four  channels  were  used  to  record  the 
voltage  on  each  cell  and  the  current  simultaneously.   This 
provided  the  experimenter  with  a  complete  picture  of  the 
behavior  of  the  battery. 

The  single  channel  chart  recorder,  a  Hewlett-Packard 
Mo s ley  680  Strip  Chart  Recorder,  Model  H01-680,  Serial  # 
712-01825  (Engineering  Projects  Lab  #1175)  was  used  to  record 
the  voltage  of  the  cell  with  the  lowest  capacity,  i.e.,  the 
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cell  which  controlled  the  end  of  discharge  point.   This  was 
cell  #3  in  battery  #2.   The  single  channel  chart  recorder  was 
used  twenty-four  hours  a  day  during  life  cycle  testing, 
primarily  because  its  slowest  speed  of  1  inch  per  hour  gave  a 
concise  record  of  charge/discharge  cycling. 

The  four  channel  chart  recorder  was  used  when  a  more 
complete  picture  of  the  charge/discharge  cycles  was  desired. 
The  slowest  speed  on  this  chart  recorder  was  5mm  per  100  sec 
(approximately  7  inches  per  hour) .   Continuous  running  would 
have  produced  an  unmanageable  amount  of  paper. 

The  set-up  for  life  cycle  testing  included  a  sensing 
circuit  to  determine  the  end  of  discharge  and  end  of  charge  and 
a  relay  to  perform  the  switching  of  the  battery  circuit,  as 
shown  in  figure  3-4.   Figure  3-5  shows  the  schematic  of  the 
sensing  and  relay-driving  circuit.   The  sensing  circuit  is  a 
Schmitt  Trigger  described  by  Millman  and  Taub  (1965) . 

It  is  a  bistable  circuit,  and  for  a  loop  gain  greater 
than  1,  its  input-versus-output  displays  hysteresis,  as 
shown  in  figure  3-6. 

The  relay  driving  circuit  simply  consists  of  a  single 
transistor  connected  in  series  with  the  relay,  with  the  out- 
put of  the  Schmitt  Trigger  to  the  base  of  transistor  Q3 .   The 
diode  connected  across  the  relay  was  simply  for  the  purpose  of 
providing  positive  opening  of  the  relay  when  the  driving 


Figure  3-5.   Schematic  Diagram  of  the  Schmitt  Trigger 

> 
Sensing  and  Relay-Driving  Circuit. 


transistor  turns  off.   The  other  relay  shown  in  the  circuit  was 
included  as  a  fail-safe  device.   It  is  energized  by  115V  ac 
line  power.   Its  terminals  are  connected  in  series  with  the 
ammeter  so  that  if  there  is  a  loss  of  power,  the  battery  cir- 
cuit is  opened.   The  terminals  of  the  main  relay  were  wired  to 
the  terminals  of  the  double  pole-double  throw  switch  through  a 
plug.   When  it  was  desired  to  perform  controlled  testing  of 
the  battery,  the  sensing  and  relay  circuit  was  unplugged. 

One  recommended  modification  to  the  circuit  is  to  insert  a 
3  position  switch  to  the  base  of  the  driving  transistor  03. 
One  position  would  connect  the  Schmitt  Trigger  to  the  tran- 
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Figure  3-6.   Hysteresis  of  Schmitt  Trigger  Circuit, 


sistor,  as  shown  in  figure  3-5.   Another  would  bias  Q3 
off,  placing  the  battery  on  discharge,  while  the  third  pos- 
ition would  turn  Q3  on  placing  the  battery  on  charge.   This 
would  facilitate  making  the  controlled  charges  and  discharges. 

tAlso,  a  simple  open  or  closed  switch  in  the  fail-safe  relay 
line  might  be  added.   With  the  set-up  as  shown ,  it  was 
necessary  to  unplug  the  relay's  energizing  power  from  the  v/all 
receptacle. 

The  Schmitt  Trigger  was  very  sensitive  to  fluctuations  of 
short  time  duration.   This  might  be  improved  by  inserting  small 
capacitors  from  the  base  of  Ql  and  Q2  to  their  emitters. 
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3.3   TESTING  PROCEDURE 

There  were  basically  two  objectives  of  the  testing. 
One  was  to  experimentally  obtain  a  family  of  discharge  curves 
for  the  verification  or  modification  of  the  Shepherd  (1965) 
equation.   The  other  was  life-cycle  testing  in  order  to 
determine  the  change  in  battery  characteristics  with  age. 

To  obtain  a  family  of  discharge  curves,  constant  current 
discharges  were  performed  at  0.6,  1.5,  3.6,  and  5.4A.   These 
values  were  chosen  because  it  was  expected  that  they  would 
provide  a  range  of  values  similar  to  those  shown  in  a  plot 
of  typical  curves,  such  as  those  in  figure  1-1.   Once  two  runs 
have  been  made  at  different  currents,  Peukert ' s  (1897) 
equation  can  be  used  to  determine  the  current  required  for  the 
discharge  times  desired. 

In  an  effort  to  place  the  battery  in  identical  conditions 
prior  to  a  discharge  in  which  the  data  would  be  used  for 
comparisons,  the  following  procedure  was  followed.   The 
battery  was  fully  charged,  followed  by  a  pre-conditioning 
discharge  at  3.6A.   The  reason  for  the  pre-conditioning  dis- 
charge was  to  remove  variations  in  the  "memory  effect"  (dis- 
cussed in  Section  1.5)  caused  by  previous  discharges.   The 
pre-conditioning  discharge  was  followed  by  a  constant  current 
charge  of  3.0A.   The  constant  current  was  maintained  by  the  use 
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Figure  3-7.   Typical  Charging  Curve  at  Constant  Current 

of  3.0A. 
of  a  12  Volt  battery  connected  in  series  with  a  variable 
resistor  and  was  used  as  a  charger. 

These  values  were  chosen  because  they  were  within  the 
range  of  currents  used  during  life  cycle  testing  and  for  con- 
venience.  By  using  this  scheme,  it  was  felt  that  repeatable 
results  could  be  obtained.   The  end  of  charge  was  determined  to 
occur  when  the  battery  was  gassing  heavily  and  the  cell  voltage 
had  remained  constant  for  about  10  to  15  minutes,  as  shown  in 
figure  3-7. 

Prior  to  discharging  the  battery,  it  was  left  on  open 
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circuit  for  a  period  of  at  least  one  hour  but  not  generally 
more  than  twelve  hours.   This  would  allow  the  electrolyte  to 
become  fairly  uniformly  distributed  throughout  the  cell  and 
also  for  the  temperature  to  stabilize. 

During  discharge,  the  variable  resistor  which  served  as 
the  load  was  adjusted  periodically  to  maintain  constant 
current.   The  end  of  discharge  occurred  when  the  voltage  of  a 
cell  dropped  off  abruptly.   Discharge  was  terminated  when  this 
voltage  dropped  to  approximately  1.0V. 

Life  cycle  testing  was  performed  using  the  Schmitt 
Trigger  circuit  and  relay  with  the  Mallory  battery  charger 
and  a  constant  load  resistance  of  about  1.6fl.   The  charging 
current  was  in  the  vicinity  of  3.0A,  ranging  from  3.5  at  the 
start  of  charge  to  about  2A  at  the  end.   The  average  discharge 
current  was  in  the  vicinity  of  3.5A.   The  Schmitt  Trigger 
circuit  was  set  to  end  discharge  and  commence  charging  when  the 
battery  voltage  dropped  to  about  5V,  and  to  switch  on  the  high 
end  at  about  7.5V.   The  setting  at  the  low  end  could  vary  by  as 
much  as  several  tenths  of  a  volt  without  any  particular  change 
in  the  time  of  switching,  due  to  the  very  rapid  fall  off  of 
voltage  at  this  point.   The  setting  at  the  high  end  needed  to 
be  much  more  precise  due  to  the  leveling  off  of  voltage  at  the 
fully  charged  point.   This  proved  to  be  most  difficult  and 
introduced  variations  in  the  state  of  charge  at  the  commence- 
ment of  discharge  under  life  cycle  testing. 
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CHAPTER  4 
RESULTS  AND  DISCUSSION 
4.1   MODIFICATIONS  TO  THE  SHEPHERD  EQUATION 


The  numeric  data  corresponding  to  figures  4-1  through  4-36 
in  this  section  is  found  in  appendix  A.   Figure  4-1  shows  the 
family  of  discharge  curves  obtained  following  the  methods  out- 
lined in  Section  3.3.   Figure  4-2  shows  the  same  curves  after 
manual  smoothing.   One  of  the  primary  effects  of  this  smoothing 
was  the  removal  of  the  initial  rapid  drop  of  voltage,  which 
Shepherd  accounted  for  in  his  activation  polarization  term,  and 
the  "coup  de  fouet."   Since  the  initial  drop  was  completed 
within  seconds  of  starting  the  discharge,  and  the  "voltage  dip" 
had  disappeared  within  about  two  minutes,  ignoring  these  pheno- 
mena should  have  little  effect  on  the  applicability  of  a  model. 

The  first  step  toward  the  establishment  of  a  model  was  an 
attempt  to  fit  the  discharge  data  to  Shepherd's  equation: 

E  -  Es  -  K  Q-^q  i    ~    Roi 

The  values  found  by  a  least  squares  fit  were: 

E  -   2.295 
s 

K  =  0.08086 

Q  =  6.844 

Ro  =  0.00092 

with  a  correlation  Z  (E  -  E    )2  for  all  curves  =  3.5008. 

expy 
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Figure    4-1.      Family   of   Experimentally   obtained   Discharge 
Curves. 
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Figure  4-2.   Smoothed  Family  of  Discharge  Curves. 
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Figure  4-3  shov/s  the  resulting  family  of  discharge 
curves.   Comparisons  of  each  individual  discharge  curve  and 
its  corresponding  modeled  discharge  curve  are  shown  in  figures 
4-4  through  4-7. 

These  results  illustrate  the  need  for  some  modification. 
The  rapid  drop-off  of  voltage  at  the  end  of  discharge  was  only 
predicted  for  the  .6  amp  discharge.   This  was  the  expected 
result,  considering  the  denominator  of  the  concentration 
polarization  term. 
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Figure  4-3.   Family  of  Curves  for  Curve  Fit  to  Shepherd's 
Equation. 
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Figure  4-4.   0.6  Amp  Discharge  Curve  and  the  Modeled  Curve  for 

the  Fit  of  the  Family  of  Curves  to  Shepherd's  Equation. 
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Figure  4-5.   1.5  Amp  Discharge  Curve  and  the  Modeled  Curve  for 
the  Fit  of  the  Family  of  Curves  to  Shepherd's  Equation. 
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Figure  4-6.   3.6  Amp  Discharge  Curve  and  the  Modeled  Curve  for 
The  Fit  of  the  Family  of  Curves  to  Shepherd's  Equation. 
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Figure  4-7.   5.4  Amp  Discharge  Curve  and  the  Modeled  Curve  for 
the  Fit  of  the  Family  of  Curves  to  Shepherd's  Equation. 


60 
Next,  a  curve  fit  for  each  individual  curve  was  per- 
formed using  the  Shepherd  equation.   Table  4-1  lists  these 
coefficients  with  the  graphical  results  shown  in  figures  4- 
through  4-11. 


Current,  i 

0.6 

1.5 

3.6 

5.4 

E 
s 

2.18 

2.18 

2.18 

2.178 

K 

0.0137 

0.0238 

0.0031 

0.00142 

Q 

6.502 

5.302 

4.373 

3.991 

Ro 

0.2779 

0.0481 

0.0734 

0.05924 

Correlation 

0.0718 

0.1239 

0.1585 

0.2638 

Table  4-1.   Curve  Fit  of  Shepherd's  Equation  to  Individual 
Discharge  Curves. 
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Figure  4-8.   0.6  Amp  Discharge  Curve  and  the  Modeled  Curve  for 
a  Fit  to  the  Shepherd  Equation  of  this  curve  alone. 
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Figure  4-9.   1.5  Amp  Discharge  Curve  and  the  Modeled  Curve  for 
a  Fit  to  the  Shepherd  Equation  of  this  Curve  Alone. 
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Figure  4-10.  3.6  Amp  Discharge  Curve  and  the  Modeled  Curve  for 
a  Fit  to  the  Shepherd  Equation  of  this  Curve  Alone. 
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Figure  4-11.   5.4  Amp  Discharge  Curve  and  the  Modeled  Curve  for 
a  Fit  to  the  Shepherd  Equation  of  this  Curve  Alone. 


63 


If  the  slope  of  the  model  were  adjusted  then  the  cor- 
relation between  discharge  and  model  would  be  very  good.   This 
suggests  modifying  the  resistance  polarization  term  such  that 
V0  =  V0(q,i)  =  (aq+b)i.   Ro  =  aq+b  takes  into  account  that  the 
ohmic  resistance  of  the  cell  changes  as  the  discharge  reaction 
proceeds  toward  completion.   This  model  simply  predicts  a 
linear  change  in  resistance.   It  may  be  found  that  R0  is  also 
a  function  of  current  since  the  PbSCU  build-up  occurs  in  a 
different  manner  at  different  discharge  rates. 

The  value  for  Q  follows  the  Peukert  equation  very  well 
when  C  and  N  have  the  values  5.803  and  1.2227  respectively. 
Figure  4-12  shows  a  log-log  plot  of  Q  versus  I  for  Peukert 's 
equation. 


Capacity 
(Ampere  -  Hours) 


Time  of 
Discharge 
(Hours ) 
100 


.1      .2     A      .6.8  1.0     2j0    k.0       6.0  &0  100 

Figure  4-12.   Log-Log  Plot  of  Peukert 's  Equation 

(i-n) 


Q  =  CI 


c  =  5.803    n  =  1.2227 


Modifying  the  Shepherd  Equation  by  the  inclusion  cf 
Peukert's  capacity  model  yielded  a  fit  which  predicted  very 
well  the  drop  in  voltage  at  the  end  of  discharge.   The 
results  obtained  did  not  produce  a  very  good  fit  to  the  left 
of  this  point,  however. 

These  plots  are  shown  in  figures  4-13  through  4-17. 
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Figure  4-13.   Family  of  Curves  for  Modification  of  Capacity 
Term. 
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Figure  4-14.   0.6  Amp  Fit  with  Capacity  Term  Modified, 
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Figure  4-15.   1.5  Amp  Fit  with  Capacity  Term  Modified. 


66 


w 
< 

En 

O 

> 


0  1 


2         3         4 


7 


Figure    4-16.       3.6   Amp   Fit  with   Capacity   Term  iYodified. 
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Figure  4-17.   5.4  Amp  Fit  with  Capacity  Term  Modified. 


67 

As  mentioned  in  Section  2.2,  the  product  of  K  and  i, 
found  from  the  curve  fits  figures  4-8  through  4-11,  tended  to 
be  constant,  as  shown  in  table  4-2. 


i         0.6  1.5  3.6  5.4 

K         0.0137       0.0238       0.00310  0.00142 

K*i       0.00822      0.0357       0.0112  0.00767 

Table  4-2.  Comparison  of  K  and  the  Product.  K*i. 


This  suggests  that  in  equation  2.5,  V-,  f   V,  (q,i),  but  rather 

V,  =V-,  (q)  =  K  -pr— - — .   By  including  both  this  new  V,  and 
da-1       Q-q  d 

capacity  terms,  the  plots  of  figures  4-18  through  4-22  were 

produced.   These  show  a  definite  improvement  over  previous  fits 

for  the  family  of  discharge  curves.  The  values  determined  were 


E   =  2.002 
s 

K  =  0.009 
R0  =  0.03 


with  a  correlation  of  1.8309. 
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Figure  4-18a.   Experimentally  Determined  Family  of  Curves 

Modification  of  V,  =  V ,  (q)  and  Capacity  are  shown  below: 
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Figure    4-18b.      Family   of   Curves    for  Modification   of 
V,   =  V , (q)    and   Capacity. 
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Figure  4-19.   0.6  Amp  Fit  with  Modification  of  V.  =  V, (q) 
and  Capacity. 
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Figure  4-20.   1.5  Amp  Fit  with  Modification  of  Vd 
and  Capacity. 
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Figure  4-21.   3.6  Amp  Fit  with  Modification  of  V,  =  V, (q) 
and  Capacity. 


2.2 

2.0 
1.8 

l.G 
1.4 
1.2 

\  1.0 

!       0.8 

0    12    3    4    5    6 

ampuH-z  kousiA 

Figure  4-22.   5.4  Amp  Fit  with  Modification  of  V,  =  V,(q) 
and  Capacity. 
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Each  plot  of  figures  4-18  through  4-22  suggests  that  if 
the  ohmic  resistance  were  modeled  as  a  linear  function  of 
accumulated  charge,  a  much  better  fit  could  be  obtained. 

Figures  4-23  through  4-26  show  the  results  of  individual 
fits  for  each  discharge  curve  by  modeling 


Ro  =  Ro  (q)  - 


Raq  +  Rt 


Table  4-3  lists  these  results 


l 

E 
s 

K 

R 
a 

Rb 
Correlation 

Figure 


0.6 

2.180 

0.00744 

0.1183 

0.03867 

0.05279 

4-23 


1.5 

2.180 

0.03453 

0.00073 

0.05424 

0.12685 

4-24 


3.6 

2.180 

0.00876 

0.01881 

0.03253 

0.03129 

4-25 


5.4 

2.180 

0.00615 

0.01771 

0.02016 

0.08875 

4-26 


Table  4-3.   Curve  Fit  for  Individual  Discharge  Curves  with 
Modifications  to  Q,  V.,  and  Ro . 
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Figure  4-23.   Fit  of  0.6  Amp  Alone  for  Modification  of 
Ro  and  Capacity. 
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Figure 
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4-24 
and 


Fit  of  1.5  Amp  Alone  for  Modif ication  of 
Capacity. 
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Figure  4-25.   Fit  of  3.6  Amp  Alone  for  Modification  of 
R0  and  Capacity. 
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Figure  4-26.   Fit  of  5.4  Amp  Alone  for  Modification  of 
Ro  and  Capacity. 
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Two  more  curve  fits  to  the  family  of  discharge  curves  were 
performed  using  the  modified  ohmic  resistance  term  and 
Peukert's  capacity  model. 

The  first  was  performed  using  Shepherd's  V,  term  with  the 
results  shown  in  figures  4-27  through  4-31.   The  second  fit 
was  made  using  the  modified  V,  term.   These  results  are  shown 
in  figures  4-32  through  4-36. 
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Figure  4-27a.   Experimentally  Determined  Family  of  Curves 

Modification  of  Capacity  and  R0  =  Ro (q)  are  shown 
below: 
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Figure  4-27b.   Family  of  Curves  with  Modification  of  Capacity 
and  R0  =  Ro (q) • 
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Figure  4-28   0.6  Amp  Fit  for  Modification  of  R0  and  Capacity, 
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Figure  4-29.   1.5  Amp  Fit  for  Modification  of  R0  and  Capacity, 
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Figure  4-30.   3.6  Amp  Fit  for  Modification  of  R0  and  Capacity. 
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Figure  4-31.   5.4  Amp  Fit  for  Modification  of  R0  and  Capacity. 
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Figure  4-32a.   Experimentally  Determined  Family  of  Curves 
Modifications  to  V-,,  R0  ,  and  Capacity  are  shown  below: 
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Figure  4-32b.   Family  of  Curves  with  Modifications  to 
V,,  R0 ,  and  Capacity. 
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Figure  4-33.   0.6  Amp  Curve  with  Modifications  -  to  V,,  Rn , 

d 

and   Capacity. 
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Figure  4-3  4.   1.5  Amp  Curve  with  Modifications  to  V,,  R0 , 
and  Capacity. 
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Figure  4-35.   3.6  Amp  Curve  with  Modifications  to  V,,  R0 , 
and  Capacity. 
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Figure    4-36.       5.4    Amp   Curve   with   Modifications    to  V^,    Ro , 
and  Capacity. 
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4.2   LIFE  CYCLE  TESTING 

As  was  previously  mentioned  in  Section  1.5,  according  to 
Bode  (1977) ,  the  ideal  requirement  for  life  testing  is  that  the 
life  test  correspond  to  the  load  conditions  of  the  battery  in 
the  field.   However,  since  ideal  testing  is  impossible  due  to 
conditions  impossible  to  duplicate  in  the  laboratory,  one 
must  choose  another  type  of  test.   Kordesch  (1977)  states  that 
the  Cycling  Life  Test  is  a  test  principally  to  evaluate  battery 
internal  components.   It  does  not  attempt  to  simulate  car 
service  conditions  which  are  extremely  variable.   Naturally, 
this  is  the  type  of  test  chosen  for  testing  the  lead-acid 
battery  which  was  used  for  experimentation  for  this  thesis. 

The  life-cycle  testing  procedure  used  to  determine  the 
change  in  battery  characteristics  with  age  is  outlined  in 
Section  3.3.   Some  of  the  pertinent  data  concerning  life- 
cycle  testing  follows: 


Number  of  Cycles  =  135 
Time  on  Charge  =19  3  hours 
Time  on  Discharge  =  114  hours 

Total  number  of  Ampere  Hours  Discharged  =  375 
Amount  of  Water  Lost  Due  to  Hydrolysis  =  12ml/cell 


To  determine  changes  in  cell  capacity  during  the  life 
of  the  battery  several  controlled  discharges  at  3.6  amps  were 
made  at  different  times  during  the  testing.   Table  4-4 
shows  these  results.   The  final  capacity  was  94%  of  the  pre- 
vious capacities.   Further  testing  is  required  to  determine 
the  rate  at  which  the  capacity  decreases . 


Number  of  Prior  Cycles   25     91     135 
Discharge  Time  72     7  2     67.5 

Amphere-hour  Capacity    4.3  2   4.32   4.05 

Table  4-4.   Capacity  Change  of  Cell  #3  During  the  Life 
of  the  Battery. 

A  great  deal  of  data  was  generated  during  the  testing  but 
it  was  considered  far  to  voluminous  to  be  included  in  its  raw 
form  in  this  report.   This  data,  which  is  primarily  in  the  form 
of  strip  chart  recordings,  was  placed  on  file  with  Professor 
H.M.  Paynter  at  the  Massachusetts  Institute  of  Technology. 

4.3   UNCERTAINTIES  IN  EXPERIMENTAL  RESULTS 

Uncertainties  involved  in  the  experimental  results  were 
due  partly  to  imperfect  knowledge  of  the  state  of  charge  and 
the  inability,  at  least  with  the  experimental  set-up  used,  to 
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control  certain  parameters.   The  most  notable  of  these  was 
temperature.   Variations  in  temperature  have  been  recognized 
to  produce  significant  changes  in  battery  characteristics. 
Great  variations  in  temperature  were  not  encountered  since  the 
experiment  was  performed  at  room  temperature,  but  there  were 
some  variations  in  this  temperature  by  nearly  as  much  as  10  C. 
The  effect  of  internal  heating  at  higher  discharge  rates  may 
have  affected  results. 

The  difficulty  of  knowing  the  precise  state  of  charge  also 
added  a  degree  of  uncertainty.   Included  in  this  was  the 
"memory  effect,"  since  the  amount  of  preconditioning  required 
to  eliminate  it  was  not  known. 

Also,  the  knowledge  of  the  mechanism  of  lead  sulfate 
build-up  at  the  various  discharge  rates  was  not  available. 

Due  to  the  differences  of  the  results  during  curve  fitting 
for  the  1.5  amp  discharge  with  respect  to  the  other  curves, 
there  was  some  suspicion  that  the  conditions  for  that  particu- 
lar discharge  may  have  been  different  than  for  the  others. 

4.4   ADDITIONAL  OBSERVATIONS 

The  method  of  charging  the  battery  and  the  point  at  which 
charging  was  terminated  was  deemed  to  be  the  most  important 
factor  in  placing  the  battery  in  the  same  state-of-charge 
before  the  controlled  discharges.   Figure  4-37  shows  a  typical 
constant  current  charging  curve  for  this  experiment. 


Voltage 
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Figure  4-37.  Typical  Constant  Current  Charge. 


The  amount  of  gassing  in  the  battery  was  determined  by 
visually  observing  the  amount  of  gas  evolution.   The  clear 
battery  case  was  vital  in  making  this  determination.   By 
observing  the  top  of  the  plates  from  the  side  of  the  battery, 
the  following  amounts  of  gassing  were  determined.   Slight 
gassing  was  determined  when  occasional  bubbles  rose  to  the 
surface.   This  was  about  the  same  as  when  the  battery  was 
standing.   Light  gassing  occurred  when  there  was  a  regular 
stream  of  these  bubbles  rising  to  the  surface.   Moderate  gas- 
sing existed  when  considerably  more  activity  was  observed. 
This  was  accompanied  on  the  charging  curve  by  a  noticeable 
increase  in  the  slope.   Heavy  gassing  was  typified  by  a  sheet 


85 

of  small  bubbles  rising  and  a  generally  more  violent  condition 
at  the  surface  of  the  electrolyte.   This  occurred  when  the 
slope  of  the  charge  curve  reached  its  maximum  and  continued 
throughout  overcharging  when  the  voltage  stopped  increasing. 

To  obtain  equivalent  charged  conditions,  the  battery  was 
overcharged  when  essentially  no  further  increase  in  voltage 
occurred  at  constant  current  for  about  15  to  30  minutes.   Most 
of  the  energy  into  the  battery  at  this  point  went  into  the 
hydrolysis  of  water,  indicating  that  the  charge  reaction  was 
essentially  completed.   The  observation  of  the  "coup  de  fouet" 
after  this  type  of  charge  and  its  absense  after  a  partial 
charge,  adds  some  validity  to  the  degree  of  completion  of  the 
reaction. 

An  explanation  of  this  phenomenon  is  that  upon  initiation 
of  discharge,  the  electrolyte  goes  into  a  meta-stable  state, 
super-saturated  with  PbSOi* ,  until    PbSCK  crystallization  sites 
start  to  form  on  the  plates.   Once  this  occurs,  the  voltage 
recovers  from  its  "dip"  and  proceeds  in  a  smooth  fashion.   This 
recovery  was  observed  to  occur  after  about  two  minutes .   The 
"coup  de  fouet"  was  not  observed  on  a  discharge  following  a 
partial  charge,  and  was  significant  for  discharges  at  higher 
discharge  rates,  3.6  and  5.4  amps,  which  followed  a  full 
charge.   The  magnitude  of  this  dip  was  approximately  0.01  volt. 

When  the  battery  was  new,  Cell  #2  had  the  lowest  capacity. 
After  several  days  of  testing,  while  moving  the  experimental 
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set-up,  Cell  #3  was  inadvertantly  placed  across  a  .10   resistor. 
After  that  occurred,  Cell  #3  had  the  lowest  capacity  and  Cell 
#1  had  the  greatest. 

After  all  testing  was  completed,  the  following  observa- 
tions of  the  battery  were  made.   The  positive  plates  and 
connectors  had  a  dark  brown  color  (artist's  color  "Burnt 
Umber") ,  which  had  a  fairly  rough  appearance.   There  was  some 
spalling  on  one  edge  of  one  of  the  plates  in  cell  #3. 

The  negative  plates  had  a  battle-ship  grey  color.   One 
small  piece  of  material  was  observed  to  have  broken  away  from 
the  plate. 

There  was  a  fairly  uniform  coating  of  sediment  in  the 
bottom  of  the  battery.   This  sludge  consisted  of  two  distinct 
materials.   One  was  a  light-colored  material  with  a  slightly 
pinkish  tint.   The  other  which  uniformly  covered  the  bottom 
of  the  battery,  and  was  layered  on  top  of  the  other  sediment, 
was  a  dark  brown  color  much  the  same  as  the  color  on  the 
positive  plates.   The  light-colored  material  may  have  been 
white  rather  than  slightly  pink  because  it  had  been  observed 
by  the  author  et  al  that  the  battery  case  itself  seemed  to  cast 
the  pink  tint  onto  the  interior  components.   It  is  believed 
that  the  light-colored  material  is  actually  predominantly 
PbS04 . 

Vinal  (1955)  states  that  overcharging  produces  sludging  of 
a  fine  brown  sediment,  and  that  prolonged  undercharging  causes 
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a  fine  white  powder  which  is  principally  PbSOi,  to  be  deposited 
on  the  bottom.   It  is  felt  that  the  dark  brown  sediment  is 
PbC>2  /  which  represents  a  loss  of  active  material  from  the 
positive  plates.   The  author  notes  that  prolonged  undercharging 
took  place  during  life-cycle  testing. 
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CHAPTER  5 
CONCLUSIONS  AND  RECOMMENDATIONS 


The  Shepherd  (1965)  equation 


E  =  Es  -  Vd(i,q)  -  V0(i)  -  Va(q)  (5.1) 

V  (q)  =  activation  polarization  =  A  exp  -£*■  (5.2) 
a                                     U 

(neglected  due  to  the  short  time  duration) 

Vo (i)  =  ohmic  polarization  =  R0i  (5.3) 

V,  (i,q)  =  diffusion  polarization  =  K  ^— f —  i  (5.4) 


proved  to  be  an  excellent  starting  point  for  establishing  an 
electrochemical  state  model  for  lead-acid  batteries.   Using  the 
lead-acid  battery  of  this  investigation,  the  author  was  able  to 
get  a  fit  to  individual  discharge  curves  which  required  only 
the  modification  of  the  ohmic  polarization  term  (equation  5.3) 
to  the  form 


V0 (i,q)  =  (Raq  +  Rfc)  i.  (5.5) 


in  order  to  get  an  excellent  fit.   The  results  obtained  from 
this  fit  suggested  that  a  better  fit  for  the  family  of  curves 
might  be  obtained  by  setting  K*i  constant. 
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The  resulting  modified  term  was 


V^  -  K  g4q  •  (5"6) 


This  was  the  form  of  the  transport  polarization  suggested  by 
Bode  (1977) .   They  also  show  that  the  use  of  the  Peukert 
(1897)  capacity  model 


Q(I)  =  CI(l_n)  (5.7) 


would  predict  the  rapid  drop  of  voltage  at  the  end  of  discharge 
very  well. 

The  modification  of  the  V-,  and  Q  terms  seemed  valid 

modeling  for  the  family  of  discharges,  but  difficulties  were 
apparent  with  the  ohmic  polarization  term.   There  were  no 

particular  consistencies  in  the  values  of  R   and  R,  for  the 

*■  a      b 

different  discharge  rates.   This  may  be  due  to  the  fact  that 
the  buildup  of  PbSOi*  at  the  plates  occurs  in  a  different  manner 
for  different  discharge  rates.   Bode  (1977)  shows  that  for 
discharge  rates  of  6  mA/cm2  and  1  mA/cm2  the  distribution  of 
PbSO^  throughout  the  positive,  capacity  limiting  plate,  is 
uniform  at  full  discharge.   For  current  densities  of  30  mA/cm2 
and  180  mA/cm2,  there  is  a  lower  PbSOi»  content  at  the  center  of 
the   plate  while  the  content  near  the  surface  is  the  same  as 
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for  the  lower  current  densities. 

The  cell  under  test  had  current  densities  of  3.4,  8.6, 
20.7,  and  31  mA/cm2  for  discharges  of  0.6,  1.5,  3.6,  and  5 . 4A 
respectively.   This  suggests  that  the  same  mechanisms  were 
encountered  in  the  battery  discharge.   Since  this  effect  is 
clearly  a  function  of  current,  a  modification  of  V0(i,q)  other 
than  that  suggested  in  equation  5.5  may  yield  a  better  fit.   It 
may  have  a  form  similar  to  Shepherd's  original  V,  term, 
equation  5.4. 

Another  consideration  might  be  that  since  equation  5.4  was 
based  on  a  small  current  density  approximation,  another  look 
might  be  made  toward  applying  equations  2.7  and  2.8  upon  which 
the  approximation  was  based. 

A  very  exciting  area  of  future  endeavor,  and  one  which 
could  prove  very  rewarding,  would  be  the  establishment  of  a 
model  which  would  predict  the  family  of  discharge  curves. 

Tremendous  advancements  in  the  lead-acid  battery  have  been 
made  since  the  time  of  Plante  and  this  progress  should  continue 
steadily  due  to  the  importance  of  this  device  in  myriad  diverse 
applications.   From  the  standpoint  of  the  System  Designer,  the 
establishment  of  a  good  state  model  is  a  very  exciting  pros- 
pect.  This  would  then  allow  for  a  determination  of  the  optimum 
use  of  the  battery,  subject  to  the  constraints  of  the  applica- 
tion. 

Further  effort  in  the  direction  of  establishing  a  model 
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which  would  predict  the  family  of  discharge  curves  should 
include  the  effect  of  discharges  at  various  rates  during  a 
single  cycle.   Other  effects  such  as  pulsed  discharge  including 
both  the  frequency  and  duty  cycle  could  then  be  added. 

The  inclusion  of  the  prediction  of  the  change  in  internal 
characteristics  during  the  secular  life  of  the  battery  is 
virtually  essential  for  an  optimization.   The  objective  in  the 
design  of  almost  any  system  is  the  optimization  of  a  cost 
(possibly  acquisition  cost  or  life  cycle  cost)  subject  to 
various  operational  criteria.   For  a  system  such  as  an  electric 
vehicle,  a  very  significant  cost  is  the  repair  cost  of  battery 
replacement.   Therefore,  it  is  necessary  to  be  able  to  predict 
the  lowest  cost  scheme  for  charge  and  discharge  battery  requir- 
ments  and  replacement  frequency  for  various  charge/  discharge 
schemes.   To  this  end,  further  life  cycle  testing  is  recommend- 
ed. 

To  produce  applicable  results,  it  is  felt  by  the  author 
that  all  cycles  should  be  identical  except  those  due  to  secular 
changes  in  the  battery.   A  rather  sophisticated  control  system 
needs  to  be  developed.   All  charges  and  discharges  should  be 
conducted  at  constant  current.   During  charge,  a  significant 
period  on  overcharge  is  recommended  to  ensure  that  the  state- 
of-charge  is  nearly  identical  prior  to  any  discharge.   Since 
the  amount  of  overcharge  cannot  be  determined  by  sensing  the 
voltage  level  alone,  some  means  needs  to  be  devised  to  deter- 
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mine  the  slope  of  the  voltage  curve  and  to  control  the  time 
on  charge  at  this  level.   It  is  also  recommended  that  a  single 
cell  be  used  during  testing  to  eliminate  any  effects  from  the 
other  cells. 

An  attempt  has  been  made  throughout  this  thesis  to  draw 
conclusions  on  the  suitability  of  the  Shepherd  (1965)  equation 
as  a  state  model  for  lead-acid  batteries.   Modifications  to 
this  equation  were  proposed,  among  them  was  the  modeling  of 
the  ohmic  polarization  as  a  linear  function  of  accumulated 
charge  to  obtain  a  good  fit  to  a  single  discharge  curve.   These 
curves  revealed  the  interesting  result  that  the  product  K*I 
was  nearly  constant.   Also,  Peukert ' s  equation  was  used  to 
model  the  cell  capacity. 

The  difficulties  encountered  in  the  attempt  to  obtain  an 
optimum  state  model  have  been  evaluated  and  enumerated.   How- 
ever/ new  discoveries  about  the  lead-acid  battery  are  being 
made,  and  new  technology  may  supply  the  means  to  overcome  the 
difficulties.   Almost  certainly,  economic,  ecological,  socio- 
logical  and  technological  problems  with  our  current  energy 
resources  will  dictate  to  society  the  necessity  of  assessing 
the  place  of  the  lead-acid  battery  in  our  future. 
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APPENDIX  A 


NUMERIC  DATA  CORRESPONDING  TO  FIGURE  k-1. 
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Numeric  data  corresponding  to  Figures  k-k    and  U-5 


Coefficients  for  the  Model 
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Numeric  data  corresponding  to  Figures  4-6  and  4-7 
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1  6  1 

0.45 

1  .  c  9  0 

1 

.  ;.  2  2 

_  r 

2  c?' 

0.90 

1  .n7  J 

1 

.7  87 

-0 . 

.  1  l-  3 

1  .3  = 

1.95  0 

1 

.74  6 

-0 . 

,2  04 

1.80 

1.920 

1 

.697 

—  !~! 

;  5  t 

C   •   £    -" 

1 .  r<  9  0 

1 

.639 

-0  . 

2  5 1 

2.7  0 

1.850 

1 

.  5  0  9 

_  0 

.  2  e  1 

3.1b 

1.800 

i 

.4  81 

—  p 

.  3 1 n 

3.4  2 

1.75Q 

1 

.417 

—  p 

■7  1  7 

3,60 

1.70  0 

\ 

.  3  6  9 

-a . 

»             L 

3.6  9 

1.660 

1 

.3  42 

-c 

.318 

3 .  7  8 

1.570 

1 

.314 

-0 

.  9  rc 

3.8  7 

1.230 

i 

.  2  8  5 

f» 

,'     ^  K 

3.9  3 

n.soo 

1 

.-c4 

~ 

,454 
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Numeric  Data  corresponding  to  Figures  k-Q    and  4-9 


Coefficients  for  the  model 


See  table  4-1 


0.6  Amp  Discharge 
Ah     Volt    Volt 
(Exp)   (Eqn) 


Diff 


1.5  Amp  Discharge 
Ah     Volt    Volt 
{Exp)    (Eqn) 


Diff 


o .  o  a 

2.110 

2.0  05 

-0. 1 05 

1  .80 

2.060 

2  .  0  0  2 

-0.C  58 

3.60 

2.010 

1  .  c95 

-  C  .  P  1  5 

4  .20 

1  .990 

1  .c*c"2 

0.000 

".80 

1  .970 

1.982 

0.012 

5.4  0 

1  .910 

1  .965 

0.0  55 

5.70 

1.880 

1  .  a  4  7 

0.0^7 

6.03 

1  .830 

1.907 

0.077 

6.12 

1.813 

1.873 

0.0  63 

fe.2* 

1.770 

1.809 

0  .  0  3a 

6.30 

1.740 

1.748 

0.008 

6.36 

1.690 

1.636 

-0 . C54 

6.39 

1.630 

1.535 

-0.C95 

6.42 

1.47  0 

1.361 

-0. 1 09 

6.4^ 

1  .030 

1.15  0 

0.120 

0  .CO 

2.052 

2  .  0  7  2 

■  ) 

'•  *"■ r. 

0.75 

2.04  0 

2 .  C  6  6 

0 

."2s 

1  .50 

2.020 

2.058 

« 

,03? 

2.25 

2.000 

?.:.*»- 

^\ 

:  it  5 

2.63 

1  . 9  ?  C 

2  .  0  3  7 

w 

,  '  4  7 

3.00 

1.980 

2 .  C  2  * 

fj 

.    as 

3.3° 

1  .963 

2.01  0 

~, 

.  "i 7 

U.13 

1.920 

1  ,?U1 

g 

.0  27 

tt.50 

1#ocr 

I.3?-1 

-0 

.03° 

4.6  5 

1.8  =  C 

1.913 

_  o 

,032 

>t.  BO 

1  .  8  1  C 

1  .73  1 

_  r* 

C79 

4. Pa 

1.77' 

1.665 

-0 

.  1  ~  5 

a .  9  c 

1.720 

1 .  c  7  0 

-•" 

.  " c  '  D 

5.0  3 

1.600 

1 .  a  2  5 

-0 

.175 

5.10 

1.10  0 

1.171 

0. 

mi 

5.13 

0  .  8  5  0 

1.03  9 

3  . 

,1  =  3 
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Numeric  data  corresponding  to  figures  U-10  and  4-11 


Coefficients  for  the  Model 


See  Table  4-1 


3-6  Amp  Discharge 
Ah    Volt    Volt 
(Exp)   (Eqn) 


Diff 


5.4  Amp  Discharge 
Ah     Volt    Volt 
(Exp)   (Eqn) 


Diff 


0 .  C  0 

2.C1-7 

1.990 

-0.127 

0  .30 

2.010 

1.389 

-0.121 

0.60 

2.0  00 

1.888 

-0.112 

0.90 

1.98a 

1.887 

-0.09-7 

1  .20 

1.970 

1*985 

-0.085 

1.50 

1  .960 

1.884 

-0.076 

1  .90 

1.945 

1.882 

-0.063 

2.10 

1.930 

1.979 

-0.051 

2.40 

1.910 

1.976 

-0.034 

2.70 

1.8  90 

1.872 

-0.019 

3.00 

1.860 

1.365 

0.00  5 

3.30 

1 . 8  3  C 

1.855 

0.025 

3.*^ 

1.790 

1.838 

0  .049 

3.73 

1.760 

1.813 

0.058 

3.96 

1  .  7  1  r 

1.793 

0.073 

a. 03 

1.55" 

1.734 

0.034 

a. 14 

1.590 

1.691 

0.101 

4.20 

1.480 

1.619 

0.139 

4.25 

1  .28° 

1.469 

0  ••  1  P.  9 

a. 32 

0.940 

0.983 

r\     nnn 

0.00. 

2.014 

1.85  C 

-0.164 

0.45 

••.900 

1.849 

-0.141 

0.30 

1.97: 

1.848' 

-0  .122 

1  .35 

1.950 

1 .  «  4  7 

-0.  103 

1.80 

1.920 

1.844 

-0.-76 

2.25 

1.890 

1.aai 

-0 .049 

2.70 

1.850 

1.834 

-  0  .  C 1 6 

3.15 

1        OA'1 

1.822 

0.072 

3  .42 

1 .75: 

1.805 

C  . '  "  5 

3.6  0 

1.70  0 

1.780 

0  .  0  °  0 

3.69 

1.66' 

1.750 

0  .  r  9  6 

3.-70 

1.57C 

1  . 7  1  3 

0.143 

3.87 

1  .23" 

1.605 

3.96 

0  .  P  r  ■" 

0.871 

0.0-71 
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Numeric  Data  corresponding  to  Figures  4-14  and  4-15 


COEFFICIENTS  FOR  THE  MODEL 


Corr 


1.872  0.00177       5.803       1.22270       3.19801 


0. 00651 


0.6  Amp  Discharge 
Ah    Volt    Volt 
(Exp)    (Eqn) 


Diff 


1.5  Amp  Discharge 
Ah     Volt    Volt 
(Exp)   (Eqn) 


Diff 


0.00 

9.110 

1.867 

-0.243 

1  .80 

2.060 

1.867 

-0.193 

3.60 

2.010 

1.866 

-0.144 

4.2: 

1.090 

1.365 

-0.125 

4  .80 

1.970 

1.364 

-0.106 

5.UC 

1.910 

1.862 

-0.0H8 

5.70 

1.88C 

1.859 

-0.021 

6.00 

1  .930 

1.854 

0.C24 

6.12 

1.310 

1 .  8  5  0 

0.040 

6.24 

1.770 

1.342 

0.072 

6.30 

1.740 

1.334 

0  .094 

5  .36 

1.690 

1.81Q 

0.  129 

6  .30 

1.630 

1.306 

0.176 

6.42 

1.470 

1.73a 

0.314 

6  .44 

1  .030 

1.757 

0.^27 

0.00 

2.052 

1 .  p  6  0 

-0. 192 

0.75 

2.040 

1.359 

-0.1=1 

1  .50 

2.02C 

1.359 

-0.161 

9    ->  ^ 

c     *    i  , 

2.000 

1  .  3  5  ; 

-0.142 

2.6  3 

1.9°'.' 

1.B5"* 

-0.133 

3  .00 

1  .gfic 

1.855 

-0  .124 

3.38 

1.953 

1.855 

-0  .  1nP 

4.13 

1.92r 

1.350 

-o.c^o 

4  .5  0 

1  .  8  8  C 

1.'^ 

-0.035 

4.6  5 

1.85: 

1.341 

-0.009 

4.80 

1.810 

1  .  3  3  '4 

0  .024 

4  .  -  o 

1  . 7  7  C 

1.829 

0  .  r-  5  3 

4.9  5 

1.72' 

1.822 

0.102 

5.03 

1.500 

1.811 

0.211 

5.10 

1.100 

1.792 

0.f^>2 

5.13 

0.8  50 

1.733 

0.933 
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Numeric  Data  corresponding  to  figures  h-l6    and  4-17 


COEFFICIENTS  FOR  THE  MODEL 


E 


K 


n 


Corr 


R, 


s  --  -o 

1.872         0.00177      5-803       1.22270       3.19808  0.00651 


3.6  Amp  Discharge 
Ah    Volt    Volt 
(Exp)   (Eqn) 


Diff 


5.4  Amp  Discharge 
Ah     Volt    Volt 
(Exp)   (Eqn] 


Diff 


0.00 

2.01"7 

1.842 

-0.17= 

0  .30 

2.010 

1.342 

-0.163 

0.60 

2.000 

1.841 

-0.159 

0.90 

1  .984 

1 .  8,  4  1 

-0.143 

1  .20 

1.970 

1.840 

-0.130 

1  .50 

1.960 

1.839 

-0.121 

1  .80 

1  .945 

1.838 

-0.1  "7 

2.10 

1.930 

1.335 

-0.094 

2.40 

1  .910 

1.834 

-0.076 

2.70 

1.890 

1.332 

-O.O 58 

3.00 

1.860 

1.828 

-0.032 

3.30 

1.830 

1.822 

-0.008 

3.60 

1.790 

1.812 

0.r22 

3.7  3 

1.760 

1.801 

0.041 

3.96 

1.710 

1.780 

0.070 

4.03 

1.650 

1.750 

0.100 

4. 1l» 

1.590 

1.724 

0.134 

4.20 

1.4  8  0 

1.578 

0.198 

4.25 

1.290 

1.578 

0.298 

4.32 

0.940 

1.19"7 

0.257 

o.oc 

2.014 

1.327 

-0.137 

0  .45 

1.990 

1 .  3  2  c 

-0.1^4 

0.90 

1.9  70 

1.525 

-0  .145 

1  .35 

1.9  50 

1.822 

-0.12° 

1  . 8  0 

1.920 

1.919 

-0.101 

2.2^ 

1.890 

1.915 

-  0  •  0  7  c 

2.70 

1.850 

1  .807 

-0.043 

3.15 

1 .  a  0 r 

1.791 

-0.009 

3.42 

1.750 

1.770 

0.020 

3.60 

1.700 

1.73  = 

0  .0  3  = 

3.69 

1.660 

1.70? 

0.048 

3.78 

1.57P 

1.652 

0  .  0  °  2 

3.37 

1  .230 

1.503 

0.279 

3.96 

0 . 8  C  C 

j. 372 

-0  .^nc> 
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Numeric  Data  corresponding  to  Figures  4-19  and  4-20 


COEFFICIENTS  FOR  THE  MODEL. 


E, 


K 


n 


Corr 


2.002  0.00900       5.803        1.22270       1.83086 


0 .03006 


0.6  Amp  Discharge 
Ah     Volt    Volt 
(Exp)   (Eqn) 


Diff 


1.5  Amp  Discharge 
Ah     Volt    Volt 
(Exp)    (Eqn; 


Diff 


0.00 

2.110 

1.975 

-0.135 

1  .30 

2.060 

1.972 

-0.0*0 

3.60 

2.01C 

1.364 

-0.046 

4.20 

1.990 

1.959 

-0.031 

a. bo 

1  .97*". 

1.950 

-0.C20 

5.  an 

1  . 9 1 C 

1.93-1 

0.021 

5.7? 

1  .  S  8  0 

1.911 

0.031 

6.00 

1.83" 

1.867 

0.037 

6.12 

1.810 

1.831 

0.021 

6.24 

1.77C 

1.761 

-0.00  9 

6.30 

1.740 

1.694 

-0.046 

6.36 

1.69C 

1.572 

-0.118 

6.39 

1.630 

1.U51 

-0.16  9 

6  .42 

1  .470 

1.270 

-0.200 

6.4  4 

1.030 

1.040 

0.010 

0 . 0  c 

1.948 

-Q  f 

104 

0.75 

2 . 0  4  C 

1.946 

-0  . 

C  "  4 

1  .50 

2.0  2" 

1.944 

-0. 

076 

2.25 

?.0Dn 

1.941 

« 

Cr° 

7.63 

1.9"" 

1  .  "  "  c 

_  f\ 

051 

3  .  C  0 

1.980 

1.936 

-0  . 

."44 

3  .3a 

1.963 

1.332 

-0  . 

.  0  3  1 

4.13 

1.9  20 

1.91* 

-0  . 

"04 

4.50 

1  .  8  »  0 

1  .ag-7 

[) 

.017 

4  .65 

1.850 

1 .  ~  9  4 

n 

,034 

4.8" 

1.810 

1.852 

0 

.0^2 

4  .83 

1.77:, 

1  .84" 

0 

.0  75 

4.9  5 

1  .720 

1  .  8  2  1 

p 

101 

5  .03 

1 . 5  :• r 

1         T9    ^ 

p 

,  1 <-. 

5.10 

1.1"  0 

1.720 

/> 

.620 

5.13 

0 . 8 c : 

1.687 

.  °  3  7 
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Numeric  data  corresponding  to  Figures  4-21  and  4-22. 


COEFFICIENTS  FOR  THE  MODEL 


K 


Corr 


2.002  0.00900       5.803       1.22270       1.83086 


0 .03006 


3.6  Amp  Discharge 
Ah    Volt   Volt 
(Exp)   (Eqn) 


Diff 


5.4  Amp  Discharge 
Ah     Volt    Volt 
(Exp)    (Eqn) 


Diff 


0.00 

2 . 0  1  7 

1.385 

-0  .132 

0  .30 

2.0V' 

1.884 

-0  .126 

0.60 

2.000 

1 .  8  8.3 

-0.117 

0.9  0 

1.9P4 

1.982 

-0.102 

1  .20 

1.970 

1.R31 

-0 .0^9 

1.50 

1.960 

1.980 

—0  .  0Q 0 

1.30 

1.945 

1  .378 

-0 .0*7 

2.1C 

1.930 

1.376 

-0.054 

2.4  0 

1.910 

1.°74 

-0.036 

2.70 

1.390 

1.8  70 

-0.02  0 

3.00 

1.860 

1.955 

0.005 

3.30 

1.830 

1.857 

0.0  27 

3.60 

1.79-1 

1.842 

0.0C2 

3. 70 

1.76  0 

1  .  5  2  ^ 

0  .06  6 

3.96 

1.710 

1.796 

0  .  0  a  6 

4.03 

1.6  50 

1.755 

0.  105 

4.14 

1.590 

1.717 

0.  127 

4  .20 

1.480 

1.55? 

0.172 

4.25 

1.280 

1.511 

0.291 

4.32 

0  .940 

0.974 

0.034 

0.0  0 

2.01a 

1 

,a31 

-0 

.18  3 

0.45 

1.900 

1 

,~3n 

—'1 

.160 

0.90 

1.970 

1 

.  9  2  8 

.142 

1  .30 

1  .950 

1 

.8  26 

-0 

.  124 

1.8: 

1  .92" 

1 

.  =  2? 

-0 

f\  r\  ■"» 

2.25 

1.890 

1 

.910 

-0 

.0-M 

2  .7  0 

1.850 

1 

.°1? 

-C 

.038 

3.10 

1.800 

1 

.7T7 

-0 

.00  3 

3.47 

1 .  7  5  r 

1 

.77£ 

r\ 

fi  ~  c 

3. 60 

1.70: 

1 

.747 

0 

.047 

3.69 

1  .66  " 

1 

.71  9 

0 

3  . 7  0 

1.57^ 

1 

.  5  6  r 

*» 

3  .87 

1.23: 

1 

.  5  3  ■'. 

'J 

3^0 

3.9  0 

0.80  ~ 

0 

.46  0 

-0 

.340 

704 
Numeric  data  corresponding  to  Figures  U-23  and  U  —  2 U . 
COEFFICIENTS  FOR  THE  MODEL 
SEE  TABLE  k-3 . 


0.6  Amp  Discharge 
Ah    Volt   Volt 
(Exp)   (Eqn) 


Diff 


1.5  Amp  Discharge 
Ah     Volt    Volt 
(Exp)   (Eqn) 


Diff 


0.Q0 

2.110 

2.102 

-0. 008 

1.30 

2.060 

2.C57 

-0.003 

3.60 

2.010 

2.009 

-0.001 

4.20 

l.Q90 

1  .991 

0.  0  01 

4.6  0 

1  .970 

1  .969 

-0.001 

5.40 

1.910 

1.940 

0.030 

5.70 

1.880 

1.916 

0.036 

6.00 

1.830 

1  .873 

0.043 

6.12 

1.810 

1  .840 

0.030 

6.24 

1.770 

1.780 

0  .010 

6.30 

1  .740 

1  .723 

-0.017 

6.3  6 

1  .690 

1  .621 

-0. 069 

6.3  5 

1.630 

1  .529 

-0.1 01 

6.4? 

1.470 

1.370 

-  C  .  1  0  0 

6.44 

1  .030 

1  .180 

0.150 

0.00 

2.052 

2.064 

Q 

.012 

0.75 

2.04  0' 

2.058 

n 

,  0  1  8 

1.5  0 

2.020 

2.049 

r 

.  0  29 

2.2C 

2.000 

2.036 

Q 

,  n  36 

2.63 

1.99  0 

2.027 

.  0  37 

3.0  0 

1  .980 

2.016 

n 

.036 

3.3  8 

1.963 

2.000 

r* 

.  0  7  7 

4.13 

1  .G20 

1  .  "  5  ? 

n 

,  r   ^  Q 

4.5  0 

1.880 

1.865 

_  p 

ft  1  ^ 

4.65 

1.850 

1.813 

_  n 

.0  37 

^.80 

1  .  8  1  Q 

1  .729 

—  p 

i  r  rtl 

4.8  8 

1.77  0 

1.664 

—  C 

]  06 

4.9  5 

1  .720 

1  .5  73 

-0  . 

.  1  <*7 

5.0  3 

1.600 

1.432 

.  1  f  8 

5  .  1  0 

1.10  0 

1.186 

0  < 

,  C"6 

5.13 

0  .850 

1.05  8 

inn 

i-      v  i  - 
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Numeric  data  corresponding  to  Figures  U-25  and  4-26 

COEFFICIENTS  FOR  THE  MODEL 
SEE  TABLE  4-3. 


3.6  Amp  Discharge 
Ah    Volt   Volt 
(Exp)   (Eqn) 


Diff 


5.4  Amp  Discharge 
Ah     Volt    Volt 
(Exp)    (Eqn) 


Diff 


0.00 

2.C17 

2.054 

0.  037 

0.3  0 

2.010 

2.033 

0  .  023 

0.60 

2.0C0 

2.C12 

0.0  12 

0.90 

1  .984 

1  .991 

0.007 

1.2.0 

1.970 

1.970 

0.000 

1.50 

1  .960 

1  .948 

-0.0  12 

1.80 

1.94  5 

1  .926 

-0.0  1° 

2.10 

1.930 

1  .904 

-0.026 

2.40 

1.910 

l.Q81 

-0. ooc 

2.7  0 

1.890 

1  .857 

-0. 0  33 

3.00 

1.860 

1.832 

-0. 028 

3.30 

1.830 

1  .804 

-0.0  26 

3.60 

1.790 

1.763 

-0.021 

3.78  ' 

1  .760 

1.741 

-0. 0 19 

3.96 

1.710 

1  .700 

-0.0  10 

<U08 

1.650 

1.652 

0.H02 

4.1<t 

1.59  0 

1  .611 

0.021 

4.20 

1.480 

1  .54* 

0.06* 

4.2  6 

1  .280 

1  .402 

0.  122 

4.3? 

0  .940 

0  .8  75 

-  0  •  0  65 

0.00 

2. 

.014 

2 

.065 

0.051 

G.45 

1 

qQn 

c 

.0  21 

r.  r  Ti 

0.90 

1 

.970 

1 

.977 

0.007 

1.35 

1 

.950 

1 

.933 

-0.017 

1.80 

1 

Q  *^>  n 

1 

.  8  8  8 

-n. 0  32 

2  .25 

i 

.890 

1 

.842 

-  0  .  fi  4  8 

2.7  0 

1 

.850 

1 

.794 

-0. 0  56 

3.15 

1 

.8  0  0 

1 

.741 

-0.0  5 Q 

3.4  2 

1 

.750 

1 

.701 

- 0 • f 49 

3  .  6  0 

1 

.703 

1 

.663 

-0. 0  37 

3.69 

1 

.660 

1 

.635 

-  0  .  n  ?  5 

3.7  8 

1 

.57  0 

1 

.591 

0.021 

3.87 

1 

.230 

1 

.  4^ 

0  •  2  6  0 

3.9  6 

Q 

.800 

0 

.750 

-  0  .  "  5  0 
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Numeric  data  corresponding  to  Figures  U-28  and  h-29 


COEFFICIENTS  FOR  THE  MODEL 


K 


Corr 


Rb 


2.020  0.00128       5-803       1.22270       3.05102       0.02100       0.00022 


0.6  Amp  Discharge 
Ah    Volt    Volt 
(Exp)   (Eqn) 


Diff 


1.5  Amp  Discharge 
Ah     Volt    Volt 
(Exp)   (Eqn; 


Diff 


0.00 

2.110 

2.019 

-0.0r-l 

1  .80 

2.060 

1  .996 

-0. 064 

3.60 

2.  010 

1.973 

-0.0  37 

4.20 

1  .99  0 

1  .965 

-0.025 

4.80 

1.970 

1  .957 

-0.0  13 

5.4Q 

1  .910 

1.943 

0  .  C  3  8 

5.70 

1.880 

1.9.42 

0.06? 

6.0  0 

1.830 

1  .935 

0.10  = 

6.12 

1.810 

1  .930 

0.120 

6.2  4 

1.770 

1.923 

0.153 

^.30 

1  .74  0 

1  .916 

0.  176 

6  •  o  b 

1.690 

1  .905 

0.2  15 

6-39 

1.63  0 

1  .895 

0.2  65 

6.42 

1  .47  0 

1  .8  73 

0.4  08 

6.44 

1.030 

1.359 

0.  829 

0.00 

2.052       7.  .  0 1 S 

-/>  .  0  I  4 

0.75 

2.040       ] 

.  =a4 

-0  .  C. 4 6 

1.50 

2  .020       : 

[  .Q70 

-0.050 

2.2  5 

2.000     : 

,  .945 

-  c . r  c  4 

2.63 

1.990        1 

.93* 

-0. 056 

3.0  0 

1.980       : 

I  .921 

U  •   U    .    . 

1.963      : 

i  .90S 

-  0  .  0  =5  £ 

4.13 

1.92  0       ] 

..881 

-0.  pt9 

4.50 

1  .880        ] 

.36  6 

U  •  i     It 

4.6C 

1.S50        ] 

.  .85^ 

0.008 

4  .  8  0 

1.810        1 

.  .  8  4  = 

0.039 

4  .  8  8 

1.77  0         1 

.  8  4  I 

0.0  73 

4.95 

1  •  7  2  0        1 

o  i  r 

0.115 

5.03 

1.600        ] 

.3  25 

0  .  2  ?5 

5.10 

1.100        1 

.3  09 

0.  709 

5.13 

0.850        1 

.301 

0.r-51 
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Numeric  Data  corresponding  to  Figures  !+-30  and  U- 31 


COEFFICIENTS  FOR  THE  MODEL 


E 


K 


Corr 


R. 


R, 


s       —  •  -      --a        "b 

2.020  0.00128       5.803       1.22270       3.05102       0.02100       0.00022 


3.6   Amp   Discharge 
Ah  Volt        Volt 

(Exp)        (Eqn) 


Diff 


5.4  Amp  Discharge 
Ah     Volt    Volt 
(Exp)    (Eqn) 


Diff 


COO 

2.017 

2.015 

-0 . 0 02 

0.30 

2.010 

1.992 

-0.018 

0.60 

2.000 

' 1.969 

-0.031 

0.90 

1  .38  4 

1  .949 

-0.038 

1.20 

1.970 

-  1  .923  ' 

-0. 0  47 

1.50 

1  .96  0 

1  .899 

-  0  .  0  6 1 

1.8  0 

1.945 

1  .876 

-0. 0  69 

2. 10 

1.93  0 

1  .852 

-0.0 78 

2.40 

1  .910 

1  .828 

-0. 062 

2.70 

1  .890 

1  .803 

-0. C  87 

3. CO 

1  .860 

1  .778 

-0. 082 

3.30 

1.830 

1  .751 

-0.  079 

3.60 

1  .790 

1.721 

-C. 0  69 

3.78 

1  .760 

1  .699 

-0.061 

3.9b 

1.710 

1  .670 

-0.040 

4.0  3 

1.65C 

1  .640 

-0.0  10 

4.14 

1.590 

1.616 

0.  026 

4.20 

1  .480 

1.578 

0.098 

4.26 

1  .280 

1  .502 

0.222 

4.32 

0  .940 

1.222 

0.2  82 

0.00 

2.014 

2 

.  G  1  2 

-  0  .  C  0  2 

0.4  5 

1  .  9C<0 

1 

.°60 

-0.070 

0.90 

1  .970 

1 

.908 

-0. 0  62 

1.35 

1.95C 

1 

.8  56 

-  0  .  P  9  4 

1.8  0 

1.92  0 

1 

.  3  0  2 

-0.3  1° 

2.25 

1.890 

1 

.7  48 

-0. 1 42 

2.70 

1.850 

1 

.6  92 

-0. 158 

3.15 

1.800 

1 

.A29 

-0.1 71 

3.42 

1.75° 

1 

.583 

-0.167 

3.60 

1  .700 

1 

.540 

-0.160 

3.69 

1.660 

1 

.508 

—  u  »  1      / 

3  .  78 

1.570 

•1 
X 

."57 

_  n  .  1 1 3 

3.8  7 

1.230 

1 

.  v-47 

C  .113 

3.9  6 

0.80  0 

J 

.510 

~  J  •  .<     -■ 

70S 

Numeric  data  corresponding  to  figures  U  —  3 3  and  i+-3^ 


COEFFICIENTS  FOR  THE  MODEL 

E        K         C       n  Corr 

s 

2.023    0.00771   5-803   .1.22270   1.391^6 

Ra  Rb 

0.015^0  0.00361 


R         R, 
a        b 

0.01540   0.00361 


0.6  Amp  Discharge 

Ah     Volt    Volt 

(Exp)  (Eqn) 


Diff 


1.5  Amp  Discharge 
Ah     Volt    Volt 
(Exp)    (Eqn; 


Diff 


2.11r 

2.013 

-r  .007 

1  .  S  0 

2.05  0 

1.994 

-0  ,ree 

3  . 6  0 

2 .  0  1  0 

tV370 

-O.C40 

4.20 

1.990 

1.9  6  0 

•0.C30 

■4.8  0 

1.970 

1.9U7 

-0.0  23 

5 .  U  0 

1  .  9 1 Q 

1  .  ?  2  c 

0.O15 

5.7  0 

1 .  arfi  c 

1.506 

0.026 

6  .CO 

1.0.3: 

1 .  S  6  5 

0.-T3S 

6.12 

1.810 

1.»33 

0.023 

5.2U 

1.770 

1.772 

0.00  2 

6.30 

1.740 

1.714 

-0.026 

6.36 

1.690 

1.509 

-0.0S1 

6.39 

1 .  6  3  C 

1.514 

-0.116 

6  .U2 

1  .47C 

1.350 

-C.120 

6.UU 

1  .  0  3  C 

1.153 

0.  123 

0 . 0  0 

1 

.05  2 

s 

."510 

-0. 

j"2 

0.75 

2 

* 

."91 

-0. 

:*;9 

1.50 

^ 

.020 

1 

.972 

-0. 

:_r 

2  .'^c; 

7 

r*  f    r\ 

1 

.9  52 

r  L.  p 

2.*3 

1 

.?r  ' 

* 

.-.40 

-0. 

r  ft  ^ 

3.0  0 

1 

.9  8C 

1 

.931 

<-■ 

:^o 

3.33 

1 

.963 

-t 

.=  1? 

-0. 

'  '15 

a.  13 

1 

.920 

1 

O  Q  O 

•                -> 

3  32 

4.50 

1 

.880 

1 

.3  63 

-0  . 

017 

a  .65 

1 

.850 

1 

.3U7 

-0. 

."  1  ^ 

4.80 

1 

.  8  1  0 

1 

.325 

0. 

015 

a .  8  8 

1 

.7  7  0 

1 

p  q  q 

C  39 

'4.9  5 

•1 
i 

.72' 

1 

.7^ 

0. 

067 

5.0  3 

1 

.  6  0  0 

1 

.754 

->  • 

154 

5.10 

1 

.10  0 

<! 

.59" 

0 . 

597 

5.13 

~; 

.  pc: 

i 

.568 

0 . 

8  1  ° 
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Numeric  Data  corresponding  to  Figures  U  —  3 5  and  U  —  36 


Coefficients  for  the  Model 

Es       K         C       n  Corr      Ra        Rb 

2.023         0.00771       5-803       1.22270       1.391^6       0.01540       0.00361 


3.6  Amp  Discharge 
Ah    Volt   Volt 
(Exp)   (Eqn) 


Diff 


5.4  Amp  Discharge 
Ah     Volt    Volt 
(Exp)    (Eqn) 


Diff 


0.0^: 

7.01"7 

''.007 

-0.015 

0  .3  0 

2.010 

1.38  5 

-0  .025 

0  .60 

2.00^ 

1.96? 

-0.032 

0.90 

1.984 

1.95r 

-0.0  34 

1  .20 

1.970 

1.93? 

-0.037 

1  .50 

1  .960 

1.915 

-0.045 

1  .80 

1.945 

1.897 

-0  .^48 

2.10 

1.93C 

1.379 

-0.r51 

2.40 

1.910 

1.360 

-0.050 

2.70 

1.890 

1.840 

-0.0^0 

3.00 

1.860 

1.819 

-0.041 

3.30 

1.830 

7.79c 

-0.035 

3  . 6  0 

1  .790 

1.766 

-0.0  24 

3.^3 

1.76  0 

1.743 

-0.017 

3.96 

1.710 

1.70"7 

-0.00  3 

4.03 

1.650 

1.56K 

0.015 

4.14 

1.590 

1.529 

0.039 

4.2C 

1.480 

1.570 

0.090 

4  .25 

1  .2°0 

1.4U* 

0.166 

4  .32 

0  . 9  4  C 

0.982 

0.042 

0  .CO 

2.014 

1 

.^" 

-0  . 

r  19 

0  .45 

1  ."9  9  0 

1 

.^5^ 

-0. 

D33 

0.9O 

1 .  9  7  C 

1 

.3  17 

-0. 

0r1 

1  .35 

1.950 

1 

.3  80 

-0. 

0  7Q 

1  . 8  0 

1  . 9  7  r, 

1 

.8  47 

-0. 

occ 

2.25 

1.890 

1 

.793 

-0. 

091 

2.70 

1  .850 

1 

.755 

-0. 

0C  5 

3.15 

1 .  9  0  r 

1 

.^0C 

n  0  c 

3.47 

1  .  7  c  0 

1 

.5  65 

-0. 

035 

3  .6° 

1.70  : 

1 

.6  2  5 

-0. 

"75 

3.6^ 

1.660 

1 

.593 

-0. 

167 

3  .  7  P 

1.570" 

1 

.540 

ft 

-  •j  •> 

3  .o"7 

1.2?" 

1 

.4  17 

0 . 
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3.96 
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